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FORTRAN PROGRAM FOR CALCULATING COOLANT FLOW AND 
METAL TEMPERATURES OF A FULL-COVERAGE- 
FILM-COOLED VANE OR BLADE 
by Peter L Meitner 

Propulsion Laboratory 
U. S. Army R&T Laboratories (AVRADCOM) 

SUMMARY 

A FORTRAN computer program called FCFC has been developed that calculates the 
coolant flow and the wall temperatures of a full -coverage -film -cooled vane or blade. 
Coolant flow is treated as one -dimensional and compressible. Heat transfer to the cool- 
ant due to impingement on the shell inner surface and convection in the film -cooling holes 
is calculated. Coolant supply pressures and main -stream gas static pressures can vary 
from hole row to row, and centrifugal effects can be included for blade calculations. 
Heat-transfer calculations can be excluded so that the program can be used as a flow 
program only. 

The vane or blade metal temperatures are calculated for the sheU inner and outer 
surfaces. All these temperatures are average values for a shell outer -surface area as- 
sociated with each fUm- cooling hole row. The heat-transfer calculations are one- 
dimensional through the wall and neglect conduction from adjacent areas. A thermal- 
barrier coating may be specified on the shell outer surface. With this option, the pro- 
gram also calculates the interface temperature between the metal and the ceramic 
coating. 

The program input is the chamber geometry (hole sizes, hole spacings, etc. ); cool- 
ant supply temperature and pressures; and main-stream gas heat -transfer coefficients, 
pressure, and velocity and temperature distributions. The physical properties of the 
coolant and the thermal conductivities of the metal and the ceramic coating are input as 
functions of temperature. The coolant flow coefficients for the impingement and film - 
cooling holes are input as functions of Mach number. The program output is a summary 
of the geometric data and the calciilated coolant -flow and heat -transfer results. 

This report presents the anal3dical procedure and identifies the necessary ass\imp- 
tions. It describes program input and output, explains error messages, illustrates two 
examples, and provides a program listing. 



INTRODUCTION 


Full -coverage film cooling is a very effective scheme for protecting turbine com- 
ponents from the hostile operating environment of high main -stream gas temperature and 
presstme. Full -coverage film cooling permits higher operating temperatures and pres- 
sures than convection cooling for greater overall cycle efficiency (lower specific fuel 
cons\imption) at acceptable coolant flow rates (ref. 1). For maximum effectiveness, 
compressor discharge air is first impinged on the inner surface of the vane or blade 
shell to remove heat by convection (ref, 2), The cooling air is then bled out through a 
large munber of evenly distribiited holes in the shell. The coolant forms a continuous, 
relatively cool, insulating layer between the shell oiiter surface and the hot main -stream 
gas. 

Numerous experiments and analyses of various aspects of full -coverage film cooling 
have appeared in the open literature, and the analysis of this report is based in part on 
the resvilts of references 3 to 6. Reference 3 approximates full -coverage film cooling as 
a form of transpiration cooling and derives the equations for metal and coolant temper- 
ature distribution for specified coolant flow, specified shell outer -surface temperature, 
and specified back -side -impingement and internal -wall heat -transfer coefficients. Ref- 
erence 4 describes a computer program that calculates the heat -transfer coefficients for 
a turbulent boundary layer on a porous wall and reference 5 describes a discrete-hole 
blowing model for full- coverage film cooling. Reference 6 establishes flow coefficients 
for a typical full- coverage-film- cooled geometry. 

Although these reports describe many aspects of full -coverage film cooling, these 
aspects have not been combined into an overall analytical procedure. Such a procedure 
has been developed and is reported herein. The coolant flow and the wall and coolant 
temperature distributions are calculated for a given vane or blade geometry; given cool- 
ant supply temperature and pressure; and given main -stream gas heat -transfer, tem- 
perature, pressure, and velocity conditions. Heat and flow balances are performed for 
each specified row of film -cooling holes and its associated portion of the shell outer sur- 
face. The flow and heat -transfer equations are solved simultaneously on the basis of 
compressible, one -dimensional fluid flow and one -dimensional heat transfer. For the 
heat -transfer calculations the equations of reference 3 are expanded and modified for a 
two-layer model to allow the inclusion of a thermal -barrier coating. Centrifugal pump- 
ing effects are included for blade calculations. 

The computer program is in FORTRAN IV and is operational on a XJNIVAC 1100/42 
computer. The program consists of 1650 cards and occupies 22 500 36 -bit words of 
memory. Its execution time is typically less than 15 seconds. 

This report explains the analytical procedvme used to develop a computer program 
called FCFC (full -coverage film cooling), which performs the described calculations. 
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The report lists the formulas used, identifies the necessary assumptions, describes the 
required program input, illustrates two examples, discusses the program outpxrt, ex- 
plains the program error messages, and provides a program listing. 


METHOD OF ANALYSIS 
Geometry and Terminology 

Figure 1 shows a section of atypical full -cover^e -film -cooled blade. Internal 
ribs, together with an insert, divide the blade cross section into individual chambers. 
The large variations in main -stream gas pressure and velocity around the airfoil periph- 
ery make chambers necessary to control and meter the coolant flow at the most advan- 
tageous local mass flux ratio, m = (pV)^(pV)g. (All symbols are defined in appendix A. ) 
The analysis described herein is for a single chamber in a vane or blade. The entire 
vane or blade is analyzed by performing the calcvilations for every chamber in that vane 
or blade. 

Figure 2 shows a cross section of a chamber and identifies the coolant flow stations. 
Station 1 is the supply plenum, station 2 is the impingement orifice plane, station 3 is 
the impingement plenum, and stations 4 and 5 are the inlet and outlet of the film -cooling 
holes, respectively. Station 6 is the main -stream gas flow immediately adjacent to the 
shell outer surface. For subsonic flow through the film -cooling holes, the static pres- 
sures at stations 5 and 6 will be equal. For sonic flow, however, the static pressime 
will be greater at station 5 than at station 6. 

The film -cooling holes in the shell are oriented by the angles a and as shown 
in figure 3. The angle a is formed by the hole centerline and its projection in the tan- 
gent plane. The angle /3 lies in the tangent plane and is measured from a chordwise 
line throTigh the hole centerline and its projection in the tangent plane. An angle of 
/3 = 0° implies in-line holes (alined in the main -stream gas flow direction), while 
j3 = 90° implies radially oriented (spanwise) holes. 


Assumptions 

The flow and heat -transfer calculations of this analysis are performed with the fol- 
lowing assumptions: 

(1) Coolant flow is one -dimensional from the supply plenum into the main -stream 

gas. 

(2) For a rotating blade, the radial pressure variation in the impingement plenum 
is that of a stationary column of air under the influence of a rotating field. 
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(3) For a rotating blade with compound -angle holes (/3 > 0 and hole entrances and 
exits at different radial locations), the press;ire changes in the film -cooling holes due to 
centrifugal pumping are much less than the normal pressure drop across the holes, 

(4) Each film -cooling hole row cools only its associated area of shell outer sxirface, 

(5) Heat transfer is one -dimensional through the vane or blade shell (stations 4 to 5). 
Calculations are performed for each specified row of film -cooling holes (includii^ back- 
side impingement and convective heat transfer in the holes), but conduction between ad- 
jacent rows is neglected. 

(6) The calculated back-side impingement heat-transfer coefficients are averaged 
over the entire inner simface (back side), of the shell. Specific impingement rows are 
not associated with specific film -cooling rows, or conversely. 


Flow Analysis 

Overall balanced coolant flow can exist throx^h a full -coverage -film -cooled chamber 
even if one or more holes have reverse flow, that is, for example, if main -stream gas 
flow travels from station 5 to station 4 in the film -cooling holes or coolant flow travels 
from station 3 to station 1 in the impingement holes. However, such a situation is un- 
acceptable from a design standpoint since any inflow of hot main-stream gas will render 
the design useless. Therefore, the flow analysis does not allow reverse flow. The de- 
tailed flow equations are presented in appendix B. 

For a given vane or blade chamber, the main -stream static pressures (station 6) 
can vary in the chordwise, as well as in the spanwise, direction. For a vane the coolant 
pressvme in the impingement plenum (station 3) will be constant, but for a rotating blade 
this pressure will vary in the radial direction. For proper flow balancing, therefore, 
each chamber in a vane or blade must be subdivided into either spanwise or chordwise 
rows of impingement and film -cooling holes, as shown in the following sections. 

Vane. - Figure 4 shows the outline of a typical full -coverage -film -cooled vane and 
one of its pressure-side chambers, which has been divided into rows of impingement 
holes and film -cooling holes with associated shell outer -siirface areas. Each shell area 
is assumed to be cooled solely by the coolant flow through the holes within that area. A 
vane impingement row consists of one or more equal -size impingement holes that have a 
common supply pressure, A vane film -cooling row consists of one or more equal -size 
film -cooling holes and the associated shell outer -surface area, which has constant main- 
stream gas temperature, pressure, and heat-transfer coefficients acting over its sur- 
face, A vane chamber can be divided into either spanwise or chordwise rows of holes, 
as illustrated in figure 4, 
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Before the coolant impingement inflow and film -cooling outflow can be calculated, 
the pressure in the impingement plenum (station 3) must be known. However, in any de- 
sign, only the supply pressures (station 1) and the main-stream gas static pressiires 
(station 6) are known. The impingement plenum pressure for balanced coolant inflow and 
outflow must be obtained in an iterative manner as follows: Avoiding reverse flow at the 
impingement and film -cooling holes requires that the impingement pleniun pressime be 
less than the lowest specified impingement supply pressure and more than the highest 
specified main -stream gas static pressure. Initially, the plenum pressvrre is taken to be 
the average of these pressTires, and the coolant inflow and outflow are calculated. If the 
resulting oxitflow is greater or less than the inflow, the plenum pressxire is decreased or 
increased, respectively, in the next flow iteration. The procedvue is continued until the 
inflow and outflow are within a relative tolerance of 0. 1 percent. 

Blade. - Figure 5 shows the outline of a typical full -coverage -film -cooled blade and 
one of its pressure -side chambers, which has been divided into chordwise rows of im- 
pingement holes and film -cooling holes with associated shell outer -siirface areas. As 
for the vane, each shell area is assumed to be cooled solely by the coolant flow through 
the holes within that area. A blade impingement row consists of one or more equal -size 
impingement holes that have a common supply pressure as well as a common radial lo- 
cation (distance from shaft centerline). A blade film -cooling row consists of one or 
, more equal -size film -cooling holes at a common radial location and the associated shell 
outer -surface area, which has constant main -stream gas temperature, pressure, and 
heat -transfer coefficients acting over its surface. A blade may be divided only into 
chordwise impingement and film -cooling rows of holes, as shown in figure 5. 

In a rotating blade, the pressures in the supply and impingement plenums (stations 1 
and 3, respectively) will vary from hub to tip. The radial supply pressure distribution 
must be specified and the resulting impingement plenum pressure distribution determined 
to calculate the coolant flow through the blade. Since there will be many rows of im- 
pingement and film -cooling holes along the span, the coolant flow from station 1 to sta- 
tion 6 will be essentially one -dimensional, with little distance traveled in the radial di- 
rection. The radial pressvire variation in the impingement plenum can thus be assmned 
to be that of a stationary column of coolant imder the influence of a rotating field. For a 
given pressure at a specific radial station (pq at station rp), the pressure at any other 
radius r is given by 


" 2/2 2\1 
-^o) 

2RTgj, 

An allowable range of base pressure is established at the minimum specified radius such 
that no reverse flow can occur at any impingement or film -cooling row. The total coolant 


p(r) = Pq exp 
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inflow and outflow are then balanced by the iterative procedure described previously for 
a vane, with the impingement plenum pressure at each impingement and film -cooling row 
calculated by the precedii^ equation. 


Heat -Transfer Analysis 

Metal and coolant temperature distributions are calculated for each shell outer - 
surface area associated with a specific film -cooling row. The detailed equations are 
presented in appendix B. These calculations cannot be done in a closed form and must 
be accomplished in an iterative manner according to the following procedure: In appen- 
dix B, the following expression for shell outer -surface temperature is obtained by con- 
sidering heat flvix through a wall: 


6 hg(0,x) - t| ihg + TjG^Cp 

This equation cannot be solved directly, since the overall effectiveness tj is also a 
function of T^ The overall effectiveness can be expressed in terms of the nondimen 
sional coolant- outlet temperature as 


V - 0c, 1^1^ - ^2 





or 






a 


2 


without and with a shell outer -stirface coating, respectively. The expressions for e - 

C j 1 

and 2 involve both the back -side impingement and film -cooling -hole heat -transfer 
coefficients. For an uncoated shell with an assumed shell outer -surface temperature 
T,,, n, the coolant temperatures at the inlet T„ . and outlet T. . of the film -cooling 

^ C^X C^O 

holes and the shell inner -surf ace temperature T - are given by 

W, 1 
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T =77(T -T )+T 
- - /V C,«> 


C,0 


Tw,l"<‘^2*‘=3)(Tw,o-Te,«> + T^,. 


For a coated shell, the coolant temperature at the film -cooling hole inlet j, at the 
interface between the metal and the coating and at the film -cooling hoie outlet 

Tg q; the shell inner -surface temperature and the interface temperature 

are given by ’ 


'^c,i ^"^w,o " "^c, 


a?' 

Co — )+Co II - — 


+ T 


C, oo 


*^c,if ^*^w,o “ "^c, 00 ^ 


Co (1 -iiV^^+co (i 


+ T 


C, °o 


T = V (T -T ) + T 
c,o '' w,o c,«>' c,«> 


"^w, i ~ ^w, o “ ^c, ^c, ®o 


Tw,« ' <^w,o - Tc. j(c2e"‘ * 

The overall iterative solution scheme is illustrated by the flow diagram of figure 6. 
Equation nximbers for cases with a thermal-barrier coating are marked with an asterisk. 
The impingement and film -cooling hole heat -transfer coefficients are fimctions of the 
calcTilated wall temperatures or coolant temperatures (thro\igh the physical properties) 
as indicated. The procedure of figure 6 is performed for every row of film -cooling holes 
at every flow-balancii^ iteration. The generated value of coolant -outlet temperature 
T^ Q affects the density and thus the calculated weight flow in the next flow iteration. 


~ PROGRAM INPUT 


The input to FCFC consists of a title card, a series of tabular input cards, and a 
series of cards describing each chamber to be analyzed. The tabular inputs are the only 
formatted data input. The data for each specific chamber are input in NAMELIST form. 
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An input data form is shown in table I, The required input cards are the title card, the 
tabular input cards, and the chamber input cards. 


Title Card 

The title card must always be present and is used to identify the particular set of 
nms. All 80 columns can be used. 


Tabular Input Cards 


The tabular input cards describe the required coolant and material physical prop- 
erties, as well as the coolant flow coefficients. Each set consists of three or more 
cards as follows: 

Card 1: NP in 12 format, where NP is the number of points in the table 
Cards 2a, 2b, 2c: the NP x-values describing the table in ascending order and in 
8F10.0 format (a maximum of 24 points) 

Cards 3a, 3b, 3c: the corresponding NP y-values in 8F10. 0 format 
The tables to be input, along with the required SI or U. S. customary units, are shown in 
table n. 

Tables 1 to 6 must always be supplied. Tables 7 and 8 can be deleted if there is no 
main -stream flow; tables 9 and 10, if no heat -transfer calculations are to take place 
(FCFC used for flow analysis only); and table 10, if there is no ceramic coating. To 
delete a table, input zero in card column 2 of the NP card. The tables of impingement- 
hole discharge coefficient (CD)., film-cooling hole total- pressure loss coefficient 
(KT)nmgj film-cooling hole flow reduction due to main-stream gas flow RT (tables 
5, 6, and 7, respectively) are given in reference 6. The program flow calculations are 
based on flow coefficients as defined in reference 6. The impingement -hole discharge 
coefficient (CD) is defined as the ratio of actual to ideal flow, the film- cooling hole total- 
pressure loss coefficient is defined as 


^^^nmg 


P3 -Ps 

Ps - Ps 


and the film -cooling hole flow reduction due to main -stream gas flow is defined as 


RT 


Actual coolant flow with main -stream gas flow 
Calculated coolant flow with no main -stream gas flow 
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The RT values of reference 6 are for a compound film -cooling hole angle ^ of 0°. 

Table 8 is used to correct RT for other values of (from 0° to 90°). 

The program FCFC generates a spline cvirve fit from each inputted set of tabvilar 
data. The curve -fitting procedure requires the slopes at the end points. These slopes 
are calculated from the first two and last two data points. For this reason, these points 
should be chosen such that fairing a straight line between them gives a good approxima- 
tion to the slope of the curve at the end points. For all tables, at least three input points 
are needed. If the program calls for a value at an x-location outside the range of the in- 
put table, the value at the nearest end point is used and an appropriate warning message 
is printed out. 

The input coordinates for table 8 are rotated through an angle of 45°, arid the spline 
fitting takes place in the rotated coordinate system. This gives a better curve fit for 
data with rapid changes in slope such as occur in input table 8. 


Chamber Input Cards 

The data for each chamber are preceded by $DATT, which is punched starting in 
card column 2. The variable names (starting in card coliunn 2 or beyond) are followed 
by an equal sign and the value or values of the variable, separated by commas. For 
each chamber, the nximber of impingement hole rows NIR and the number of film -cooling 
hole rows NFCR are specified; the maximvtm allowable rows are 25 and 50, respectively. 
Subscripted variables are associated with specific rows; that is, the subscripted 
value is associated with the row of holes. When fewer than the maximum number of 
rows are specified, subscripted variables need only have as many input values as the 
specified nximber of rows. Integer values must be input without decimal points. The 
last data value for each chamber is followed by a $ instead of a comma. The input data 
are retained for multiple chamber inputs. Thus, if a variable is common to successive 
chambers, it has to be input just once for the initial chamber. The chamber geometry 
input variables are defined by figime 7. All chamber input variables, along with the re- 
quired SI or U. S. customary vmits, are shown in table HI. 

The variables lUNTS to OMG in table m specify the types of calculations desired. 
These variables have been assigned default values as shown. The variables NIR to RGAS 
are associated with the impingement hole rows: NIR is the number of specified impinge- 
ment hole rows, and NIHPR to PIT are subscripted variables associated with the im- 
pingement rows. As such, each variable must have at least NIR input values. The var- 
iable HSPl, the hole spacing for each impingement row, is used in determining the back- 
side impingement heat -transfer coefficient (eq. (BID). This correlation is based on a 
square impingement array, with equal spacings in the spanwise and chordwise directions, 
as shown in figure 7. In practice, however, these spacings may differ and the average 
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of the two spacings should then be specified. The variables TT and EGAS define the 
coolant gas; they are not subscripted and are thus constant for all rows of impingement 
holes. 

The variables NFCR to ROV2G of table HI are associated with the film -cooling hole 
rows. The variable NFCR is the number of specified film -cooling hole rows, and 
NFCHPR to ROV2G are subscripted variables that must have at least NFCR values. The 
variable HSP5 is the hole spacing for each film -cooling hole row (fig, 7), and,' as for 
HSPl, an unequal array spacing should be reduced to an equivalent square spacing. The 
variable HFC4 (h factor at station 4; fig. 2) is a modification factor for the calcxilated 
impingement heat -transfer coefficient at each film-cooling hole row. For the film- 
cooling heat -transfer calculations, the calculated impingement heat -transfer coefficients 
are averaged over the shell back side (inner sm-face), since the program does not asso- 
ciate specific impingement rows with certain film -cooling -hole rows, or conversely. 
When back -side heat -transfer coefficients vary (from centrifugal effects or from im- 
pinging at less than perpendicular to the surface), HFC4 is a multiplier used to modify 
the back -side heat -transfer coefficient at the specified film -cooling -hole rows, (This 
variable has a default value of 1. 0. ) The variable HFC45 (h factor for stations 4 to 5; 
fig. 2) is a mviltiplier used to modify the calculated film-cooling hole heat-transfer co- 
efficients for each row (eq. (B13)). Equation (B13) is valid for hole length -diameter 
ratios L/D between 1.0 and 8.0. For L/D less than 1.0, reference 7 measured heat- 
transfer coefficients that were as much as 50 percent greater than predicted by equa- 
tion (B13) (entrance effects). The correction factor HFC45, which has a default value of 
1. 0, is used to account for this. The variable TMSG is the main -stream gas temper - 
atxure, which must be the same as the temperature used to evaluate the main -stream gas 
heat -transfer coefficients. 


PROGRAM OUTPUT 

The FCFC output is a printout of the title card, the input data for all specified tables, 
and the calculated results for each chamber. The chamber output consists of the follow- 
ing messages and blocks of tabxilated data; 

OUTPUT FOR CHAMBER XX 


Units and Option Messages 


XX ROWS OF IMPINGEMENT HOLES 


Impingement Hole Input Data 
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XX ROWS OF FILM COOLING HOLES 


Film- Cooling Hole Input Data 

IMPINGEMENT AND FILM COOLING FLOWS HAVE CONVERGED IN 
XX OVERALL ITERATIONS 

INFLOW EQUALS XXXXX.XXX KG/HR (LBM/HR) 

Impingement Flow Results 

OUTFLOW EQUALS XXXXX.XXX KG/HR (LBM/HR) 


Film- Cooling Flow Results 
HEAT TRANSFER RESULTS 
Heat- Transfer Results 

Each of these blocks is described in the following subsections. 


Units and Option Messages 

One or more of the following messages about the system of units and the particular 
options used are printed out; 

SI (ENGLISH) SYSTEM OF UNITS 

COOLANT GAS CONSTANT = XXXXXX.XX J/(KG-K) ((FT-LBF)/(LBM-R)) 

THIS CASE IS FLOW ANALYSIS ONLY AND INCLUDES NO METAL TEMPERATURE 
CALCULATIONS 

THIS CASE INCLUDES A THERMAL BARRIER COATING 

THIS CASE INCLUDES CENTRIFUGAL EFFECTS. ROTATIONAL SPEED ECUALS 
XXXXX.XX RPM 
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Impingement Hole Input Data 



Film -Cooling Hole Input Data 

This block of output tabulates the following for each row of film -cooling holes: 
ROW row number 

HOLES number of holes per row 

DIAMETER hole diameter, mm; in. 

THICKNESS 

WALL wall metal thickness, mm; in. 

COATING coating thickness, mm; in. 

L/D hole ler^h -diameter ratio 

HOLE SPACING hole spacing, mm; in. 

ALPHA hole chordwise inclination angle 

BETA hole compound inclination angle 

2 2 

RHOVG main -stream gas value of pV, kg/m "hr; Ibm/ft -hr 

2 2 2 

RHOV2G main -stream gas value of pV , kg/m- hr ; Ibm/ft-hr 
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R4 distance from shaft centerline, mm; in, 

2 

P6 main -stream gas static back pressure, N/cm ; psia 

The L/D is that value associated with the combined thickness of the wall and any speci- 
fied coating. When no main-stream flow is specified (MSBL=0), main-stream gas 
RHOVG and RHOV2G are printed as zero. The variable R4 is the location of the film- 
cooling hole centerline on the shell inner surface. For noncentrifugal calculations, R4 
is printed as zero. 


Impingement Flow Results 

This block of output tabulates the following for each row of impingement holes: 


IMP ROW 

row number 

PSPLYT 

2 

coolant supply total pressime, N/cm ; psia 

P2 

2 

static pressure, N/cm ; psia 

M2 

Mach number 

T2T 

total temperature, K; °F 

T2 

static temperature, K; °F 

WIMP 

coolant inflow, kg/hr; Ibm/hr 

CDIMP 

impingement discharge coefficient 


The coolant supply total pressure, shown as PIT in the section Impingement Hole Input 
Data, is repeated here as PSPLYT, 

Film -Cooling Flow Results 

This block of output tabulates the following for each row of film -cooling holes: 


FC ROW 

row number 

P3T 

2 

impingement plenum pressvtre, N/cm ; psia 

P4 

2 

static pressvtre at inlet, N/cm ; psia 

M4 

Mach number at inlet 

T4T 

total temperattme at inlet, K; °F 

T4 

static temperature at inlet, K; °F 

P5T 

total pressure at exit, N/cm ; psia 
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P5 

M5 

T5T 

T5 

TCTIF 

WOUT 

KT 

RT 

RT 

CORR 

RHOV 

RATIO 

RHOVSQ 

RATIO 

ITRS 


2 

static pressure at exit, N/cm ; psia 
Mach number at exit 
total temperature at exit, K; 
static temperature at exit, K; °F 

total coolant temperature at metal -coating interface, K; °F 

coolant outflow, kg/hr; Ibm/hr 

total-pressure loss coefficient 

reduction in coolant flow due to main -stream flow 

correction factor for RT 

ratio of coolant -to -main -stream density times velocity 

ratio of coolant-to -main -stream density times velocity squared 

number of iterations needed to achieve film -cooling flow convergence 
in last overall flow iteration 


When no coating is specified (KCLC=0), the coolant interface total temperature prints 
zeros. When no main-stream flow is specified (MSBL=0), the pV and pV ratios print 
zeros and RT and RT CORR print 1, 0. The main-stream pressure, shown as P6 in the 
section Film -Cooling Hole Input Data, is repeated here as P5. K the flow through the 
film -cooling holes is subsonic, P5 and P6 will be equal. However, for choked flow, P5 
will be that pressure determined from the compressible -flow relations at Mach 1. 0 and 
will be greater than the specified main -stream pressure P6. 


Heat -Transfer Results 

This block of output tabulates the following for each row of film -cooling holes: 

FC ROW row number 

HEAT TRANSFER 
COEFFICIENTS: 

HGO main-stream gas heat -transfer coefficient for coolant temperature 

2 2 

equal to main-stream gas temperature, J/m -sec • K; Btu/ft ■ 
hr -°R 
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HGl 


main -stream gas heat -transfer coefficient for coolant temperature 

2 2 

equal to shell outer-surface temperature, J/m • sec • K; Btu/ft • 
hr - °R 

2 

FC-HOLE heat -transfer coefficient in film -cooling hole, J/m • sec-K; 

Btu/ft^ . hr • °R 

2 

IMPG back -side impingement heat -transfer coefficient, J/m -sec-K; 

Btu/ft^ . hr • °R 

H MODIFICATION 
FACTORS: 

FC-HOLE modification factor for film -cooling hole heat -transfer coefficient 

(inputted HFC45) 

IMPG modification factor for back -side impingement heat -transfer coeffi- 

cient (inputted HFC4) 

2 2 

COOLED AREA cooled area associated with each film -cooling row. cm ; in. 

GAS TEMP main-stream gas temperature, K; *^F 

WALL TEMP- 
ERATURE: 

OUTSIDE shell outer -surface temperature, K; °F 

INTFACE shell interface temperature, K; °F 

INSIDE shell inner -surface temperature, K; °F 

AVG. THERM. 

COND. : 

METAL metal average thermal conductivity, J/cni- sec- K; Btu/ft- hr- °R 

COATING coating average thermal conductivity, J cm- sec-K; Btu/ft- hr- °R 

ETA overall effectiveness 

ITR number of iterations required to achieve metal temperature con- 

vergence in last flow iteration 

The tabulated values of film -cooling hole and impingement heat-transfer coefficients 
include their corresponding modification factors. When no coating is specified (KCLC=0), 
the interface temperatures and coatir^ thermal conductivities are set to zero. The aver- 
age thermal conductivities for the metal and coating are evaluated at the average temper- 
atures through the metal and coating, respectively. 
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Error Messages 


Error messages have been incorporated in the calculation procedure. The messages 
for the main program and the various subroutines, along with possible causes and cor- 
rective actions, are as follows (where error messages that do not stop program execu- 
tion are preceded by the word "WARNING'’): 

Main program . - The error messages for the main program are 

CASE ABORTED - A REQUIRED CURVE WAS NOT INPUT OR WAS SPECIFIED 
BY LESS THAN 3 POINTS 

Check the required input tables and add the missing data or specify at least three points, 
CASE ABORTED - COATING WAS SPECIFIED BUT NO COATING THICKNESS 
Specify coating thickness. 

CASE ABORTED - THE SPECIFIED PRESSURES WILL RESULT IN REVERSE 
FLOW 

Check the specified supply and back pressures or alter hole sizes, 

WARNING - T2 HAS NOT CONVERGED IN 15 ITERATIONS FOR IMPINGEMENT 
ROW XX 

This message could be caused by specifying significantly erroneous physical properties. 

WARNING - T5 HAS NOT CONVERGED IN 15 ITERATIONS FOR FILM COOLING 
ROW XX 

WARNING - T5T HAS NOT CONVERGED IN 15 ITERATIONS IN OVERALL 
FLOW ITERATION XX 

These messages could be caused by specifying significantly erroneous physical prop- 
erties or the heat -transfer -coefficient modification factor HFC45. 

WARNING - THE AVERAGE PRESSURE BETWEEN STATIONS 4 AND 5 HAS 
NOT CONVERGED IN 15 ITERATIONS FOR FILM COOLING ROW 
XX 

WARNING - P5T HAS NOT CONVERGED IN 15 ITERATIONS FOR FILM 
COOLING ROW XX 

These messages could be caused by specifying significantly erroneous physical prop- 
erties or the total-pressure loss coefficient curve 

IMPINGEMENT AND FILM COOLING FLOWS HAVE NOT CONVERGED IN 
25 ITERATIONS 

Change hole sizes and/or supply and back pressures. 
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Subroutine TMETO. - The error message for subroutine TMETO is 


WARNING - OUTER WALL TEMPERATURE HAS NOT CONVERGED IN 15 
ITERATIONS IN OVERALL FLOW ITERATION XX 

This error message can be caused by specifyii^ erroneous values of the main-stream 
gas heat -transfer coefficients HGO and HGl. The message can also be caused by the 
initial values assumed in the iterative process. If the message appears for values of 
overall flow iteration that are less than the actual number of flow iterations required 
(given by the message "IMPINGEMENT AND FILM COOLING FLOWS HAVE CON- 
VERGED IN XX OVERALL ITERATIONS"), the solution is valid. 

Subroutine MNEW . - The error message for subroutine MNEW is 

WARNING - M HAS NOT CONVERGED IN 25 ITERATIONS 

Check the inputted table of total- pressure loss coefficient (KT) . 

nrng 

Subroutine SPLINE . - The error messages for subroutine SPLINE are 

WARNING - A SPECIFIED X -VALUE (XXXXX.XXX) IS BELOW THE RANGE 
OF INPUT TABLE XX 

WARNING - A SPECIFIED X-VALUE (XXXXX.XXX) IS ABOVE THE RANGE 
OF INPUT TABLE XX 

Check the inputted tables and extend their range as required. 


EXAMPLE PROBLEMS 

The use of FCFC is illustrated by analyzing a chamber on both the vane and blade of 
a high -temperature, high-pressure core turbine. Example 1 demonstrates flow and heat- 
transfer calculations for a vane chamber with a thermal -barrier coating. Example 2 
demonstrates centrifugal flow calculations without heat transfer and thus shows how 
FCFC can be used as a flow program only. 

The inputted tables of impingement discharge coefficient (CD)., film -cooling hole 
total-pressure loss coefficient (KT)j^g, and film -cooling hole flow reduction due to 
main-stream gas flow RT were obtained from reference 6. The main-stream gas heat- 
transfer coefficients HGO and HGl were evaluated by using the Stanford University 
STAN5 computer program of reference 4 which was modified to include the discrete-hole 
blowing model of reference 5. 
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Example 1 


A section of the vane and chamber that were analyzed is shown in figure 8. Also 

shown are the impingement hole diameters; the film -cooling hole diameters; the main- 

2 

stream gas pressures P6; and the associated main -stream gas values of pV, pV , HGO, 
and HGl. The vane material is MAR-M509 and the coating is yttria -stabilized zirconia 
(Y202-Zr02). The vane span is 3. 81 centimeters (1. 50 in. ), and the impingement and 
film -cooling hole spacings are 0. 381 and 0. 254 centimeter (0. 15 and 0. 10 in. ), respec- 
tively. The shell and thermal -barrier coating thicknesses are 0. 127 and 0. 0127 centi- 
meter (0. 050 and 0. 005 in. ), respectively, with an impingement distance of 0. 0889 centi- 
meter (0. 035 in. ). Coolant supply pressure is 404 N/cm (586 psia) and coolant tem- 
perature is 811 K (1000° F). Main-stream gas hot-spot temperature is 2550 K 
(4130° F). 


Example 2 

A section of the blade and chamber that were analyzed is shown in figure 9. The 
blade span and the impingement and film -cooling hole spacings are the same as for the 
vane of example 1. Impingement and film -cooling hole sizes are constant at 0.4318 and 
0.4572 millimeter (0.017 and 0.018 in.), respectively. For this example, which in- 
volves no heat -transfer calculations, the wall thickness and the impingement distances 
were taken to be constant at 1. 016 and 0. 762 millimeter (0. 040 and 0. 030 in. ), respec- 
tively. In the actual blade, both vary from hub to tip. Coolant supply temperature was 
811 K (1000° F). The analysis was further simplified by assuming an impingement and 
film -cooling row at each of 15 specified radial locations. (In general, impingement and 
film -cooling rows are staggered. ) Also, each film -cooling row was taken to consist of 
two adjacent holes (one from each chordwise station) and was assumed to have a radial 
position equal to the average radial position of the two holes (fig. 9). The radial varia- 
tions of coolant supply pressure PIT and main-stream gas values of static pressure P6, 
2 

pV, and pV for the 15 rows are tabulated in figure 9. 

Table IV lists the input data for the two example problems. The title card, the tab- 
ular inputs, and the chamber inputs are identified. Tables V to VII show the program 
output for the two examples. Table V shows the title card and all tabular data. Tables 
VI and Vn are the outputs for the vane and blade chambers, respectively. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 22, 1978, 

505-04. 
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APPENDIX A 


A 

CD 

S 

^1 ■ ^6 
D 

F 

G 

Sc 

h 

KT 

k 

L 

I 

M 

m 

N 

Nu 

Pr 

P 

q 

R 

Re 


SYMBOLS 

area, ft^ 

parameters defined by eqs. (C5) and (C6) 
discharge coefficient 

specific heat at constant pressure, J/(g-K); Btu/(lbm-°R) 
constants of integration defined by eq. (C45) 
diameter, m; ft 

function values at specified input points 

flow rate per unit area, kg/(m •!«:); lbm/(ft -hr) 

force-mass conversion constant, 1; 32. 174 (lbm)(ft)/(lbf)(sec ) 

porous -wall -matrix, internal, volumetric, heat -transfer coefficient defined 
by eq. (CIO), J/(m^-hr.K); Btu/(ft^- hr- °R) 

heat -transfer coefficient, J/(m^*hr*K); Btu/(ft^* hr- °R) 

porous -wall -matrix, heat -transfer coefficient, J/(m *hr-K); Btu/(ft -hr- R) 

total- pressure loss coefficient for flow into still air 

thermal conductivity, J/(m-hr-K); Btu/(ft* hr- °R). 

thickness, m; ft 

length, m; ft 

Mach number 

blowing ratio, (pV)^/(pV)g 

dimensionless heat -transfer -coefficient parameter defined by eq. (C17) 

Nusselt number 

Prandtl number 

pressxire, N/m^; Ibf/ft^ 

heat fl\ix, J/(m^-hr); Btu/(ft^-hr) 

gas constant, J/(kg- K); ft-lbf/(lbm-°R) 

Reynolds number 
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RT ratio of coolant flow with main-stream gas flow to coolant flow without main- 
stream gas flow 

r radius, m; ft 

T temperature, K; °R 

V velocity, m/sec; ft/sec 

W flow rate, kg/hr; Ibm/hr 

X distance, m; ft 

y function value at any arbitrary ordinate location 

Z porous -wall -matrix, internal sirnface area per unit voliime, 1/m; 1 /ft 

a film -cooling hole inclination angle, deg 

a^, Og coefficients defined by eqs. (C30) and (C31) 

^ film -cooling hole compound angle; or parameter defined by eq. (C9) 

y ratio of specific heats 

77 overall effectiveness, defined by eq. (B16) 

0 dimensionless temperature parameter defined by eq. (B18) 

\ parameter defined by eq. (C 8 ) 

II kinematic viscosity, kg/(m-sec); lbm/(ft-sec) 

^ dimensionless distance parameter defined by eq. (C7) 

p density, kg/m^; Ibm/ft^ 

<p dimensionless temperatvme parameter defined by eq. (B17) 

parameter defined by eq. (C44) 
u) rotational speed, 1 /sec 

Subscripts: 

a based on impingement -jet arrival velocity 

av average 

b bulk 

c coolant 

ct coating 

fc film cooling 

g main -stream gas 
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i inner surface 

if interface 

imp impingement 

loc local 

m metal 

n based on impingement hole centers 

nmg no main-stream gas 

o outer surface 

w wall 

0 base 

1 station at supply plenum 

2 station at impingement orifice 

3 station at impingement plenum 

4 station at film -cooling hole inlet 

5 station at film -cooling hole exit 

6 station at shell ovrter surface in main -stream gas flow 

oo free stream; or supply 

Superscript: 

’ total 
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APPENDIX B 


EQUATIONS 
Flow Equations 

Impingement flow. - The coolant flow rate through the impingement holes (treated as 
orifice flow) is given by 

Wimp = <CD)impP2V2Ai„p (Bl) 

where 




Film cooling flow. - The coolant flow rate through the film -cooling holes (treated as 
pipe flow with friction) is given by 

Wjp = (B4) 

where 
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p. _ P3 + Psg^Wg 

1.0 + (KT)„^g 


(B 7 ) 


Mach number change across a film -cooling hole. - Consider a constant film -cooling 
hole area. When the hole exit Mach number and the total temperature and pressure, as 
well as the change in total temperature and pressure across the hole are known, the hole 
entrance Mach number can be obtained as follows: If the inlet station is designated by 
subscript 4 and the outlet station by subscript 5 , the continuity equation gives 


P4V4A4 - P5V5A5 


(B8) 


Equation (B8) can be expressed as 


p^M^A^ ^/y^RT; 




rt; (1 + 


^4 


- 1 


(74+1)72(74 -1) 


M‘ 


RTg 1 1 + 




- 1 


(75+1)72(75-1) 


M, 


(B 9 ) 


Solving for M4 gives 


rt; 1 1 + 


r/1 - 1 


(74+1)72(74-1) 


•M" 


M 4 = 


If 1^5 ~ ^ 2 

RTg fl + -2 M^ 


(75+1)72(75-1) 

5 1 P4A4 


(BIO) 


This equation is solved iteratively by Newton's method. 


Heat -Transfer Equations 

Back-side impingement . - The heat-transfer coefficient on the shell inner surface is 
calculated from the Garden -Cobonpue impingement correlation (ref. 8) 


0 . 286 k(Re)?-® 2 ^ 


av 


(BID 
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Convection in film -cooling holes. - The heat -transfer coefficient in the film -cooling 
holes is calculated from the Davey correlation (ref. 7), from which the local Nusselt 
number varies along the length of the hole as 


-0 2/T 

(Nu)ioc = 0. 036(Re)®- ^ j j 


(B12) 


From the definition of Nusselt number, the average heat -transfer coefficient over the 
entire length of the hole I is obtained by integrating 


/ 


I 


hoc 


av 


= 0 


/T 0 2 

. 045 - (Re)°- ®(Pr)°- ( —\ * 

D \t^) \lj 


(B13) 


The average heat -transfer coefficient in the portion of the hole between stations I ^ and 


1 2 is evaluated from 


/ 


*^av = 




hjQ^dx 0.045(^)(Rer-“(Pr) 


^0.8,^.0A 'b 


0. 18 


D 


0.2 


w/ 




I2 -li 




(B14) 


Shell outer -surface temperature. - Heat flux through a wall can be expressed as 


q = n) = G C„(T, n - T. «,) = GCV (T^ , - T^ 


g' g w,o' “c p' c,o ^c, 


c p ' ' w,o c, 


(B15) 


The overall effectiveness rj is defined by 


V 


T - T 
_ c,o c, °° 

T - T 
w,o c, °° 


(B16) 


After we introduce the parameters 


<P 


T - T 
^ g w,o 

T - T 

g c,°° 


(B17 
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and 


e = 



-T 

-T 


c,o 

w,o 


equation (B17) can be reduced to 


(p = 


hg + GpCp’’ 


By assvuning constant properties and using superposition (ref, 9), 

hg(e,x) = hg(0,x) - e[hg(0,x) - hg(l,x)] 


(B18) 


(B19) 


(B20) 


or 


hg(0,x) = hg(0,x) - G Ahg (B21) 

where hg(0,x) and hg(l,x) are the heat -transfer coefficients for the coolant temperature 
equal to the gas temperature and the shell outer -surface temperature, respectively. 
These heat -transfer coefficients are obtained from a suitable boundary -layer computer 
program and are based on an initially assximed shell outer -surface temperature. 

The dimensionless temperature groxipings can be combined to give 

e = r ~ (B22) 

<p 

Combining equations (B19), (B21), and (B22) then gives 


hg(0,x) - ?? Ahg + 7GpCp 


(B23) 


This equation can be solved for T^ ^ to give 


T 

w,o 



(Tg-T^,j[pG^Cp-f(l-p)AhJ 

hg( 0 ,x) -77 Ahg + T/Gj.Cp 


(B24) 


y 


25 



Full -coverage film cooling . - Consider the cross section of a coated, full -cover age - 
film -cooled wall as shown in sketch (a). 



The coolant temperatui'es are designated by T at the supply, T . at the film- 

c , ^ ^ 

cooling hole inlet, T .. at the interface between the metal and the coating, and T 

Co j 11 C/ ^ vJ 

at the film -cooling hole outlet. The metal temperatures are designated by T . at the 

W, 1 

shell inner surface, T^ at the interface between the wall and the coating, and 
T at the shell outer surface. The main -stream gas temperature T_ is that temper- 

W j O g 

ature in terms of which the main -stream gas heat -transfer coefficients are evaluated. 
Reference 3 develops an analytical model to predict the coolant temperature rise and 
the metal temperatvtre distribution through a porous wall. The results hold for fixed 
values of shell outer -sxurface temperahme coolant temperattire T^ and im- 

pingement and film -cooling hole heat -transfer coefficients. For a single metal layer, 
the coefficients resulting from the specified boimdary conditions can be solved for ex- 
plicitly. The solution takes the form 


where 


and 


=Ci +036 +Cge ^ 


(B25) 


0^(0 +C 2 U +Co fl 


(B26) 




T - T 
w c, °° 

T - T 
w,o c,°° 


(B27) 


(»c(5) 


T - T 
c c, °° 

T - T 
w,o c,°° 


(B28) 
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are the nondimensionalized temperature distributions in the wall and coolant, respec- 
tively. All symbols are defined in appendix C where the analytical model for a two -layer 
wall is also developed. The equations for each layer take the same form, but the six re- 
sulting constants cannot be solved for explicitly and must be evaluated numerically. The 
solution is 




(B29) 


«c, 1(^1^ = - C2 ( 1 - . C3 (1 - 


(B30) 


V2(^2)=^4^C5e ^ ^ ^ C^e ^ 2 


(B31) 


2 \ / 2 ^ 

®l \ / Ot^ \ 0^rt ^ rt 

^c,2(^2) = C4-.C5^1-^le ^Kc,il-A]e 


(B32) 


The subscripts 1 and 2 on 0 and 0 refer to the metal and coating, respectively. The 

W 

constants Cj, C 2 , and Cg for the two -layer wall are different from the corresponding 
constants for the one -layer wall. 

The overall effectiveness 77 is given by 


and 



1 ( 1 ) - C 2 





(B33) 


7j - 0^ 2 ( 1 ) ~ ^4 + ^ t 




(B34) 


for an uncoated and a coated shell, respectively. For an uncoated shell, T^ ^ T^ 

and T are given by 
w, o 



w,o 


- T ) 
c, °o' 





+ C 3 



(B35) 
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(B36) 


Tw,i = ^"^2+C3)(T^^o-Tc,»^^^c,oo 

For a shell with a thermal -barrier coating, j, •, and 

evaluated from 


T . = (T - T ) 
c,i ' w,o c,“° 


4' 

Co |l - — 1+Co 11 --i 


+ T 


C, °o 




'^ 2 !' --I" 


a« \ a2 


^ly 


+ T 


c.°° 


^C,0 ^W,0 " ”^C,' 


’'w,1 = '^w,o-Tc.<,KC2*C3) + T^_„ 


Tw,lI = <Vo-’^c,«Hv''-C3e''KT 


C, <=o 


(B37) 

are 

(B38) 

(B39) 

(B40) 

(B41) 

(B42) 
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APPENDIX C 


DERIVATION OF EQUATIONS FOR METAL TEMPERATURE DISTRIBUTION AND 
COOLANT TEMPERATURE RISE IN A TWO-LAYER POROUS WALL 

Reference 3 develops the equations for metal temperature distribution and coolant 
temperature rise through a single -layer porous wall with a fixed shell outer -surf ace tem- 
perature. The results are 


and 


where 






+ Cge 




(Cl) 


ee(^)=Ci + C2 




(C2) 


0 = Z 

W,0 “ c,«> 

(C3) 

„ 

0 =— 2 ^ 

V rp _ T’ 

w,o " c,<» 

(C4) 


(C5) 


(C6) 

^ =- 
L 

(C7) 

H L^ 
X- “ 

(C8) 


k 
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(C9) 


/3 = 


m 

G„C„ 
c p 


= hZ 
m m 


The boundary conditions are shown to be 


eji) = 1 


N0^(O) = e;(0) 


0c(o) = f ^;(0) 

A 


and the constants of integration are 


Cj = 0 


N 


^2 ■ 




a< 

(N - a2>e ^ - (N - aj^)e 




where 


^3 - 


aj - N 


(N - a2)e - (N - aj)e 


N = 


h.L 


Now consider a two -layer porous wall as shown in sketch (b). 


(CIO) 

(Cll) 

(C12) 

(C13) 

(C14) 

(C15) 

(C16) 


(C17) 
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Let the shell outer -surface temperature be 
nate the inner and outer layer, respectively. 


Q and let the subscripts 1 and 2 desig- 
Using the equations 


^i=- 


so - — 


iLi 

X. = 

•L 1.. 


h ~ 


oLo 
m,z i 


/3 = ^ 

^ G C 
c p 




(C18) 

(C19) 

(C20) 

(C21) 

(C22) 

(C23) 


results in the following wall temperature and coolant temperahire expressions for each 
layer; 


V I<h> = Cl + 


(C24) 


«c, 1«1> = Cl * Cj * Cj 


(C25) 


and 


®w,2^^2^ ■ ^4 ^5® ^R® 


“ 1^2 . ^ “ 2^2 
6® 


(C26) 
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e 


(C27) 


where 




"2^2 



\ + 4 ^ 1 ] 

(C28) 


)j-y/(i2+4Xj) 

(C29) 


^2 "^^2 “^^2 ) 

(C30) 



(C31) 


The six constants are evaluated from the boundary conditions as follows: As in ref- 
erence 3, an energy balance at boundary 1 leads to 

Ni0w,l(O) = 6w,l(0) (C32) 

and 

0c i(0) = — 0^ l(0) (C33) 

^1 

At the interface between the two layers (boundary 2) there must be continuity in metal and 
coolant temperatures, as well as continuity in heat flux. This is expressed by 


and 



(C34) 

9c, 1<1> = 9c, 2<®' 

(C35) 

9 ;, 1 <«=^ 9 ;_ 2 ( 0 ) 

(C36) 
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Finally, at boundary 3, the specified wall temperature gives 
Substituting equations (C24) to (C27) into equations (C32) to (C37) then gives 


NjCj + (Nj^ - aj^)Cg + (N^ - 32 ^^3 ~ ® 


C1 + C2 



C j + Cg© 


+ 036^2 = C4 + C5 


+ C, 


Cl + Cg 





= C4+C5 



(C38) 

(C39) 

(C40) 
^ (C41) 


where 


^1 ^2 

naje Cg + C3 - OjCg - OgCg = 0 


Q!. QU 

C4 + Cge ^ + Cge ^ = 1 


n = 


4^1 


(C42) 

(C43) 


(C44) 


From equation (C39) it can be shown that Ci = 0, Other than that, no further simplifi- 
cation is possible and the remaining constants (Cg to Cg) are best solved by a matrix 
solution from 


33 





(C45) 
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APPENDIX D 


PROGRAM STRUCTURE AND FUNCTION 

The computer program FCFC consists of the main program MAINP and the sub- 
routines TMETO, MNEW, AIRPRP, PRBMTX, SPLINE, and XMTXSL. The calling 
relations between MAINP and the subroutines are shown in figure 10. The functions of 
MAINP and each of the subroutines are described in this appendix. 


Main Program MAINP 

The main program MAINP is the control program that directs the flow of the solution 
from input to output and calculates and balances the coolant flow. Program MAINP reads 
the input, makes the necessary conversions to working units, establishes the initial 
plenum pressure or pressure profile (for centrifugal calculations), balances the coolant 
outflow and inflow by an iterative procedure, prints the output, and returns the variables 
to the input units. Flow and heat transfer are solved simultaneously, with all heat- 
transfer results being obtained from the TMETO subroutine. 


Subroutine TMETO 

Subroutine TMETO performs all heat -transfer calculations including back -side im- 
pingement, convection in the film -cooling holes, and full -coverage film cooling. It cal- 
cvilates the heat picked up by the coolant at all flow stations and the inner and outer tern - 
peratures of the metal and the thermal -barrier coating. 


Subroutine MNEW 

Subroutine MNEW establishes the Mach niimber at the inlet of a constant -area film- 
cooling hole, for a given total temperature and pressure at the hole exit, and for a given 
change in total temperature and pressure across the hole (eq. (BIO)). 


Subroutine AIRPRP 

Subroutine AIRPRP calculates the physical properties of the coolant at any specified 
temperature. The properties are evaluated from input tables 1 to 4 by calling subroutine 
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SPLINE. Subroutine AIRPRP performs any necessary imit conversions (from SI into 
U. S. customary imitsj and calculates values of different combinations of gamma: y - 1, 
(y - l)/y, y + 1, (y + l)/2, y/(y - 1), and (y - l)/2. The Prandtl number is evaluated 
from its definition Pr = C pA- 


Subroutine PRBMTX 

Subroutine PRBMTX evaluates the function second derivatives at the specified x- 
locations for all input tables. The slopes at the end points are evaluated from the first 
two and last two data points. The calculation of the second derivatives was separated 
from the spline -fitting procedure of subroutine SPLINE, since the second derivatives 
have to be calculated only once but the spline -fitting procedttre is performed many times. 


Subroutine SPLINE 

Subroutine SPLINE generates an interpolated (spline fitted) value of y at any x for 
a curve described by a finite number of points (ref. 10). 


Subroutine XMTXSL 

Subroutine XMTXSL is a general matrix -solution technique based on the Gauss 
Jordan elimination method (ref. 11). 
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APPENDIX E 


PROGRAM VARIABLES DICTIONARY 

The variables used in the main program and in the subroutines are described here. 
Subscripted variables pertaining to the impingement and film -coolir^ rows are shown 
with the indexes I and J, respectively. Variables that are input arguments in a sub- 
routine are defined in the listing of the calling program. 


A5(J) 

AIMP(I) 

ALPHA(J) 

ANEW 

ANGRl,... 

ANGRIO 

AOS 

AOLD 

AOUT(J) 

BETA(J) 

CDI(I) 

CDD 

CDFC(J) 

CDIFC 

CDOD 

CFFLOW 

CFMCH 

CFP45 

CFP5T 


Main Program MAINP 

shell outer -surface area associated with the film -cooling row 

impingement -row hole area 

film -cooling -row inclination angle 

hole area at entrance of film -cooling hole (dummy variable for constant - 
area hole) 

rotation angle for coordinate system of input tables 1 to 10 
input argument for AOUT(J) in subroutine TMETO 

hole area at exit of film -cooling hole (dummy variable for constant -area 
hole) 

film -cooling -row hole area 

film -cooling -row compound angle 

impingement -hole discharge coefficient 

output argument for curve 5 in SPLINE subroutine 

film -cooling flow reduction due to main-stream blowing 

temporary storage for CDI(I) 

output argument for CDFC(J) in subroutine SPLINE 
relative tolerance for total inflow and outflow 

relative tolerance for Mach number iteration between stations 4 and 5 
relative tolerance for P45 
relative tolerance for P5T 
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CFT2 

CFT5 

CFT5T 

CFTWO 

CP 

DAU 

DAU2 

DFC(J) 

DI(I) 

FCBLR 

FCHD 

FCHSP 

FLOFC 

G 

GAM 

GCVG 

GDGMl 

GMl 

GM1D2 

GMIDG 

GPl 

GP1D2 

GTST 

HO 

HI 

HFC4(J) 

HFC45(J) 

HFCTR 

HG0(J) 


relative tolerance for T2 
relative tolerance for T5 
relative tolerance for T5T 

relative tolerance for shell outer -surface temperature 

specific heat at constant pressure 

input argument for TAU{J) in subroutine TMETO 

input argument for TAUC(J) in subroutine TMETO 

film -cooling -row hole diameter 

impingement -row hole diameter 

film -cooling blowing rate (input argument in subroutine TMETO) 

input argument for DFC(J) in subroutine TMETO 

input argument for HSP5(J) in subroutine TMETO 

relative chaise between total coolant inflow and outflow 

specific -heat ratio, r 

y evaluated at next -to -last value of TN 

relative change between GTST and GAM 

y/(y - 1) 

y - 1 

(r - D/2 

(y - l)/y 
y + 1 
(y + l)/2 

y evaluated at last value of TN 

input argument for HGO(J) in subroutine TMETO 

input argument for HGl(J) in subroutine TMETO 

modification factor for impingement h 

modification factor for film -cooling hole convective h 

input argument for HFC4(J) in subroutine TMETO 

main -stream gas h for coolant temperature equal to main -stream 
gas temperature 
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HGl(J) 

HHFCTR 

HSP 

HSPlfn 

HSP5(J) 

ICTR 

IHLD 

IJ 

lOA 

lUNTS 

JCV(J) 

JCVT 

JHLD 

JRVFL 

JRVFLT 

K 

KCLC 

KCNVG(J) 

KKLM(J) 

MSBL 

MTC 

NC 

NFCHPR(J) 

NFCR 

NIHPR(I) 

NIR 

NPCl,... , 
NPCIO 

NREAD 


main -stream gas h for coolant temperature equal to shell outer -surface 
temperature 

input argument for HFC45(J) in subroutine TMETO 

ratio of film -cooling hole spacing to diameter 

impingement hole spacing 

film -cooling hole spacing 

indicator for centrifugal calculations 

indicator for supply row with lowest specified R1 

coimter for overall flow iterations 

coxmter for chamber calculations 

indicator for SI or U. S. customary units 

convergence indicator 

chamber convergence indicator 

indicator for film -cooling row with lowest specified R4 

film -cooling reverse -flow indicator for individual rows 

film -cooling reverse -flow indicator for entire chamber 

counter for overall film- cooling flow iterations 

indicator for coating or no coating 

counter for individual film -cooling flow iterations 

counter for individual film -cooling -row heat -transfer calculations 

indicator for main -stream gas blowing 

indicator for metal temperature calculations 

input table number 

number of film -cooling holes per row 
number of film -cooling rows 
number of impingement holes per row 
nxunber of impingement rows 

number of points specified for input tables 1 to 10 
integer number of input read file 
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NWRITE 

OMG 

PIT (I) 

PITHLD 

PITMIN 

P2(I) 

P2T(I) 

P3T 

P3TFK 

P3TFCR(J) 

P3TIR(I) 

P3TMNTSr 

P3TMNR 

P3TMXX 

P4(J) 

P4T(J) 

P45 

P45CNV 

P45HLD 

P45N 

P45T 

P5(J) 

P5HOLD 

P5MAX 

P5T(J) 

P5TCV(J) 

P5TNEW 

P5TOLD 


integer number of output write file 

rotative speed 

total pressure at station 1 

temporary storage location for PIT 

minimum specified supply pressure 

static pressure at station 2 

total pressure at station 2 

total pressure at station 3 (vane calculations) 

temporary storage for P3T 

total pressure in impingement plenum at each film -cooling row 
(blade calculations) 

total pressure in impingement plenum at each impingement row 
(blade calculations) 

lowest allowable pressure in impingement plenum 

total pressure in impingement plenum at minimum specified radius 

highest allowable pressure in impingement plenum 

static pressure at station 4 

total pressure at station 4 

average static pressure in film -cooling hole 

relative change in P45 

next -to -last iterated value of P45 

last iterated value of P45 

average total pressure in film -cooling hole 

static pressure at station 5 

temporary storage for P5 

highest specified back pressure for vane calculations 

total pressure at station 5 

relative change in P5T 

last iterated value of P5T 

next -to -last iterated value of P5T 
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P6(J) static pressure at station 6 

PFCR temporary storage location for P3TFCR(J) 

PHOLD temporary storage location for P3TMXX or P3TMNN 

PN45(J) static pressure at midpoint of film -cooling hole 
PRN Prandtl number 

PTN input argument for P4T(J) in subroutine MNEW 

PTO input argument for P5T(J) in subroutine MNEW 

Rl(l) radial distance at station 1 

R4(J) radial distance at station 4 

REJ2(I) Reynolds number at station 2 

REJ5(J) Reynolds number at station 5 

REYN45 Reynolds number at midpoint of film -cooling hole 
RGAS gas constant 

RIHLD temporary storage location for R1(I) 

R4HLD temporary storage location for R4(J) 

RH02(I) density at station 2 

RH04(J) density at station 4 

RH045 density at midpoint of film -cooling hole 

RH05(J) density at station 5 

RMN lowest specified R1(I) or R4(J) 

RIMN lowest specified R1(I) 

R4MN lowest specified R4(J) 

ROV2C(J) pV^ of coolant at station 5 
ROVG(J) pV of main -stream gas 

ROV2G(J) pV^ of main -stream gas 

ROV2R input argument for ROV2RT(J) in subroutine SPLINE 
ROVRAT(J) (pV)^/(pV)g 
ROV2RT(J) (pV^) /(pV^L 

RTCOR output argument for RTCR(J) in subroutine SPLINE 
RTCR(J) correction factor for CDFC(J) 
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T2(I) 

T4(J) 

T45 

T5(J) 

TAU(J) 

TAUC(J) 

TAUI(I) 

TC 

TC2(I) 

T2CNVG 

T5CNVG 

TD 

T2D 

T5D 

TERM 

TG 

T2HLD 

T5HLD 

TITLE 

TMI(J) 

TMO(J) 

TMSG(J) 

TN 

TO 

TT 

T2T(I) 

T4T(J) 

T5T(J) 

T3TAV 

T4TAV 


static temperature at station 2 

static temperature at station 4 

average static temperature between stations 4 and 5 

static temperature at station 5 

shell metal thickness 

shell coating thickness 

impingement insert thickness 

input argument for TT in subroutine TMETO 

coolant interface temperature (boundary 2) 

relative change for T2(I) 

relative change for T5(J) 

temporary storage for T4(J) or T4T(J) 

input argument for T2(I) in subroutine AIRPRP 

input argument for T5(J) in subroutine AIRPRP 

|l.O - [P2a)]/P2T(I)j^'>'“l^/'>' or |l.O - [P5(I)]/P5T(I)y^"^^/’' 

input argument for TMSG(J) in subroutine TMETO 

temporary storage location for T2 (I) 

temporary storage for T5(J) 

title of calculations 

inner -wall temperature 

outer -wall temperature 

main -stream gas temperature 

output argiunent for T4(J) in subroutine MNEW 

input argument for T4(J) in subroutine MNEW 

coolant total supply temperature 

total temperature at station 2 

total temperature at station 4 

total temperature at station 5 

average coolant total temperature at station 3 

average coolant total temperature at station 4 
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TTN 

TTO 

T5TFTR 

T5TOLD(J) 

TW2(J) 

V2(I) 

V4(J) 

V45 

V5(J) 

WFCR 

WIMP (I) 

WIMPT 

WOUT(J) 

WOUTT 

XBETA 

XCDI 

XDI 

XETA(J) 

XHD(J) 

XHH(J) 

XHSPl 

XILOD 

XIMP(I) 

XKA 

XKT 

XKTD 

XLC 

XLFC(J) 

XLFCC(J) 

XLFCPC(J) 


input argument for T4T(J) in subroutine MNEW 

input argument for T5T(J) in subroutine MNEW 

relative change in T5T(J) 

next -to -last iterated value of T5T(J) 

wall interface temperature 

velocity at station 2 

velocity at station 4 

average velocity in film -cooling row 

velocity at station 5 

input argument for WOUT(J) in subroutine TMETO 

coolant inflow 

total coolant inflow 

coolant outflow 

total coolant outflow 

input argmnent for BETA(J) in subroutine SPLINE 

average impingement discharge coefficient 

average impingement hole diameter 

overall effectiveness 

impingement heat -transfer coefficient 

heat -transfer coefficient in film -cooling holes 

average impingement hole spacing 

ratio of impingement distance to impingement hole diameter 

impingement distance 

coolant thermal conductivity 

film -cooling total -pressure loss coefficient 

output argument for curve 6 in SPLINE s\jbroutine 

inp;it argiunent for XLFCC(J) in subroutine TMETO 

length of film -cooling hole (metal only) 

length of film -cooling hole (coating only) 

length of film -cooling hole (metal plus coating) 



XLM 

XLODFC(J) 

XLODI(I) 

XLODXX(J) 

XM2(I) 

XM4(J) 

XM5(J) 

XMD 

XMK1(24),... 

XMK10(24) 

XMNEW 

XMOLD 

XMU 

XRH02 

XT4TAV 

XV2 

XXI,... , 
XXIO 

XXAKCT(J) 

XXAKM(J) 

XXIMP 

XXKT(J) 

YYl,..., 

YYIO 

ZFC 


A1 

A2 

AKCT 


input argument for XLFC(J) in subroutine TMETO 

film -cooling hole length -diameter ratio (metal only) 

impingement hole length -diameter ratio 

film -cooling hole ler^h -diameter ratio (metal plus coating) 

Mach number at station 2 

Mach nvimber at station 4 

Mach number at station 5 

temporary storage location for XM2(I) or XM5(J) 

’ calculated values of cvirve slopes Mj^ for tables 1 to 10 

output argument for subroutine MNEW 

inpvit argument for XM5(J) in subroutine MNEW 

coolant viscosity 

average density at station 2 

average total temperature at station 4 

average velocity at station 2 

X -coordinates for input tables 1 to 10 

coating thermal conductivity 
metal thermal conductivity 
average impingement distance 
total -pressure loss coefficient 

y -coordinates for input tables 1 to 10 

input argument for XLODFC(J) in subroutine TMETO 

Subroutine TMETO 

parameter defined by eq. (C5) 
parameter defined by eq. (C6) 
coating thermal conductivity 
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AKM 


metal thermal conductivity 


ALl 

AL2 

AREAR 

BETA 

BETA2 

C2,... ,C6 

CMAT(24,25) 

CN(24) 

COEF 

DA 

DA2 

DELHG 

DEN 

ETA 

FACVA 

HC 

HD 

HH 

HH2 

HM 

HM2 

KLM 

REH 

RENA 

ROOT 


parameter defined by eq. (C30) 
parameter defined by eq. (C31) 
area reduction ratio 
parameter defined by eq. (C22) 
parameter defined by eq. (C23) 
constants obtained by solving eq. (C45) 

general problem matrix to be solved by subroutine XMTSOL 
solution vector obtained from subroutine XMTSOL 
coefficient (temporary storage location) 
parameter defined by eq. (C20) 
parameter defined by eq. (C21) 

HO - HI 

denominator of eqs. (CIS) and (C16) 

overall effectiveness, defined by eqs. (B33) or (B34) 

arrival velocity factor 

HD corrected for presence of film -cooling holes 

coolant impingement -heat -transfer coefficient obtained from Gardon- 
Cobonpue correlation (eq. (Bll), ref. 8) 

average convective -heat -transfer coefficient in film -cooling hole 
(metal only, eq. (B13)) 

average convective -heat -transfer coefficient in film -cooling hole 
(coating only, eq. (B14)) 

internal volumetric -heat -transfer coefficient (metal only) 
internal volumetric -heat -transfer coefficient (coating only) 
counter for number of wall temperature calculation iterations 
Reynolds number in film -cooling hole 

impingement Reynolds number based on "arrival" velocity 
+ 4X^ 


ROOT2 



TCA 

TCAO 

TCCAV 

TCIF 

TCIN 

TCO 

TCTAV 

TDIF 

TFILM 

TNEW 

TOLD 

TR 

TWAV 

TWI 

TWIF 

TWO 

TWOCVG 

U 

XLTOT 


CNVCR 

DNM 

I 

PATG 

POWN 

POWO 

XMFCN 

XMFCO 


average coolant temperature in film -cooling hole (metal only) 
overall average coolant temperature 

average coolant temperature in film -cooling hole (coating only) 

coolant temperature at interface plane 

coolant temperature at inlet of film -cooling hole 

coolant temperature at outlet of film -cooling hole 

coating average temperature 

temperature difference, TG - TC 

/ 

film temperature, (TWI + TC)/2 
last iterated value of TWO 
next -to -last iterated value of TWO 
temperature ratio 

average wall temperature (metal only) 

wall inner temperature 

wall interface temperature 

wall outer temperature 

relative change in TWO 

parameter defined by eq. (C17) 

total length of film -cooling hole (metal and coating) 

Subroutine MNEW 

relative tolerance for Mach number iteration 
denominator in expression for iterated Mach number 
counter for Mach number convergence iteration 
(P5/P4KA5/A4) 

(y + l)/[2(y - 1)] evaluated at last value of y 

(y + l)/[2(y - 1)] evaluated at next-to-last value of y 

1. 0 + [(y - 1)/2]M^ evaluated at last value of M 

O 

1. 0 + [(y - 1)/2]M evaluated at next-to-last value of M 
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XMHLD 

XMN 

XMNEW 

XNUM 


temporary storage location for XMN 
Mach number 

final iterated value of XMN 

numerator in expression for iterated Mach number 


ANGROT 

CAN 

F(24) 

FR(24) 

L{24) 

LR(24) 

MAT(24,25) 

N 

NPl 

NP2 

OPT 

SAN 

SOL (2 4) 

XK(24) 

XKR(24) 

XPFST 

YPFSTR 

YPLST 

YPLSTR 


Subroutine PRBMTX 

coordinate system rotation angle 
cos (ANGROT) 

specified points that describe curve in unrotated coordinate system 
generated points that describe curve in rotated coordinate system 
lengths of intervals between inputted F(24) in unrotated coordinate system 
lengths of intervals between generated FR(24) in rotated coordinate system 
matrix of function second derivatives at specified XK locations 
number of intervals generated by XK values of FR (NPl - ]) 
number of points that describe a curve. N + 1 
N + 2 

indicator for rotated coordinate system 
sin (ANGROT) 

solution vector of problem matrix MAT (24,25) 

inputted x -values corresponding to inputted points F(24) in unrotated 
coordinate system 

generated x -values for a rotated coordinate system 
slope of first interval in unrotated coordinate system 
slope of first interval in rotated coordinate system 
slope of last interval in unrotated coordinate system 
slope of last interval in rotated coordinate system 
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Subroutine SPLINE 


ANGLNV 

ANGROT 

CAN 

CANT 

CRIT 

DELXR 

FK 

FKMl 

FR(24) 

IND 

LK 

MK 

MKMl 

N 

NC 

NMl 

OPT 

SAN 

SANI 

TERMl,... , 
TERM4 

X 

XKR(24) 

XR 

XX 

XXM 

Y 


inverse of coordinate system rotation angle (-ANGROT) 
coordinate system rotation angle 
cos (ANGROT) 
cos (ANGINV) 

relative accuracy of iterated y -value for rotated coordinate system 
(XX - XXM)/10 

value of specified function at first point to right of desired x-location 

value of specified function at first point to left of desired x-location 

specified y -values of table in rotated coordinate system 

indicator for determining whether desired x-value is outside inputted 
range of x 

length of interval 

value of function second derivative on right side of interval 
value of function second derivative on left side of interval 
number of intervals that describe a curve 
input table number 
N - 1 

indicator for rotated or unrotated coordinate system 
sin (ANGROT) 
sin (ANGINV) 

terms whose sum re^equal to spline -fitted value of y 

x-location in unrotated coordinate system 

specified table x -locations in rotated coordinate system 

x-location in rotated coordinate system 

specified x-value on right side of interval 

specified x-value on left side of interval 

spline -fitted value at specified x in unrotated coordinate system 
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Subroutine XMTXSL 


DET 

DIV 

FCT(24) 

ISNGL 

MAT (24, 

NC 

NLST 

NM 

NN 

NR 

NSW 

SOL(24) 
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matrix determinant 
value of row pivot element 

factor used to reduce elements in pivot column to zero 
factor for indicating singular matrix 

overall matrix obtained by adding problem matrix and identity matrix 
number of columns 
NC +NR 
NR - 1 
NC + 1 

number of rows (order of matrix) 

number of switches needed to make pivot element the largest element 
solution vector 
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APPENDIX F 


PROGRAM LISTING 


H(1N PR06RAH 


OINCNSiaN TITLEII6I 

DIMENSION NIHPRI25I,Rl(?SI,0II25»,T»Ur(?5),HSPl(2S),XIMP(?5» ,PIT(? 
»5> 

DIMENSION NECHPRISOI.RmSOt.OrCCSOItASiSOl.TAUtS^I.HSPSrSIl.HFUMrS 
«0),HFCNS tSOI.ALPH* ISO! .BETAISOI ,HEQ t SD > *HG1 1 SO > ,T MSE I 50 > . »b ( SD ) ,R0 

• VE( 5O>,R0VZG(50t ,TAUC(S0l,TU2IS0l,TC2f SO > 

DIMENSION AIMPC2S) , XLOOI « 25 1 ,P3TI R f ?5 1 , XM2 1 25 ) , X2 1 25 » , 12 ( ?5 » , T2 T ( ? 

• 5 I,P2(25I ,P2H25I,C0I (25 I ,RH02 1 25» , RE J2 • 25 1 ,U IMP ( 2 5 1 

DIMENSION AOUTISOI ,*LFC ( SO I .XLOOFC ( 50 » ,P 3TFC R (50 > , JC V ( 53 1 ,KCNVG(50 
*),XM5(50»,NS«50J,T5(5OI,T5Tf50l,PS(5O»,P5T(5O»,T5T0L0(50I,C0FC(5D) 
*,XXXT(S0»,RHO5(5QI ,ROVRAT(50I,ROX2C(5O),ROi(2RI(50l ,RE J5 ( 5 0 » , XL F CC ( 
•501 .XLFCPC(SQI.XLOOXX (50t.RICR(50l 
DIMENSION TAISOI.mTISOl.PNISOl.PilKSOl.VlKSOI.RNIlStSOl.IMKSOl.TM 

• 0 (50 I ,XETA(501,XMA (SO I ,UOUT( 50 1 .P5TC V ( 50 1 ,RHO4(50 I 
DIMENSION XHD(5a>,XHH(50»«XXAKM(501,KKLM(50)tXXAKCT(501 
DIMENSION XXI (24 I.XX2(2NI ,XX3I2<I>,XX4(2N 1,XX5(24I ,XX5(2in .XX7(2n I , 

»XX8 (241 ,XX9(241 ,XX 10(24 I 

DIMENSION Yyi(241,YY2(24l,Yr3(24),YV4(24l,YY5(24»,TTb(24»,*Y7l24), 
*YY8(241,YY9(241,YYI0(24I 

DIMENSION XMK1(241,XMK2(24 1 ,XMK3(24 ),XMX4(24),XMX5(24),X1<b(24l,XM 
*K7( 24 ),XHK8 (24 > ,XNK9( 24 I ,XHKtO(24 ) 

COMMON NPCl .NPC2,NPC3tNPC4,ANSRI.«NER2t ANGR3t*N5R4 ,XX1 ,XX?,XX3,XX4 
• ,YY 1 ,YY2,YY3,YY4 ,XMX1 , XMK2 , XMK 3 t XMK4 ,NRE AD , N WRI T E 

NAMELIST/OATT/1UNTS,1CTR,MTC,MSBL,XCLC,OMG,RGAS, 

•NIR,N1HPR,R1,0I ,TAUI,HSP»,XIMP,P1T,TT, 

*NFCR,NFCHPR,R4,0FC,A5.TAU,TAUC,HSP5,HFC4 ,HFC4 5 , AL PHA , 5 E I A ,HGO , HG 1 , 
*TMSG,Pb,RO¥G,ROV2G 

NREAD=5 

NKRlIE^b 

READ (NRE A 0,50 101 (TITLEdltl-l.Ibl 
HPITE (NURITE.bOlO) (TI TLE ( I I , I- 1 , 1 b I 

ANGR1=0. 

READ(NREA0,5020IN»C1 
IFCNPCl .LT. 31GOTO 130 
RFAO(NREA0,5033 I (XXI ( I ) ,Irl ,NPC1> 

READ(NREA0,5030) (YYH I ) ,Irl ,NPCI 1 
UPITE(N¥RITE, 60201 
DO ID I-1,NPC1 

13 WR11E(NWRITE,6030»(XX1 (Il.YYllI 11 

CALL PH8MTX(NPC1,XX1,YYI,ANGR1,XMX1 1 

ANGR2rO. 

READ(NREAD,5020)NPC2 
1F(NPC2 .LT. 3IGOTO 130 



RE»0IMRE»D,50JOI (*X2« !>•*=* 

RC<0(NRE«6,5030l(rY?Cl I ,l=l,NPC?» 

UR1TE{NHRITE,60R0I 

DO 20 I:l.NPC2 

20 URITEfNURlTE,6QS0IIXX2(II»yV2CIII 
C»Lt PRBRTXINPC2,XX2, VY2,»N6R2,XRK2I 

ANGR3=0. 

Rr«0(NREk0«502DINPC3 

IFINPC3 .LT.3IG0T0 130 

READ (NREA0,5Q30 MXX3( 1 1,1=1. NPC3I 

READCNREA0,5030> (rY3(II,I=l.NPC3i 

URITEINURITE, 60601 

00 30 I=I ,NPC3 

30 UPITEINWR1TE,6030I (XX3(II .YY3IIII 
CALL PRBHTX INPC3,XX3. YY3.ANGR3.XRN3 I 

ANGR<I=0. 

REAO(NREAO,‘>02D)NPC<I 
IMNPCR .LT. 3IG0T0 130 
READCNREAD,S030> (XXRt I t ,I=1,NPC4I 
REAO(NREAD,S030I (YY<I( I > ,I = l,NPC>i> 
URITEINURITE, 6070> 

00 PO 1=1 ,NPCP 

AO UPITE(NURITE.60S0I (XXPdl.YYPdtl 

CALL PRBNTXINPCR.XXP.YYP.ANGRP.XNKPI 

ANGR5=D. 

RCAO(NREAD,5020INPCS 
IMNPCS .LT. 31G0T0 130 
RCA0(NREA0,S030) ( X XS ( 1 t ,1= 1 ,NPC5 I 
READ tNRE A 0,5030 I (YY5( I> ,1 = ] .NPC5I 
UPITE (NUR1TE,60BO) 

00 50 1=1 ,NPC5 

50 UPITE(NURITE,6030,I <XX5(I» ,YY5dd 
CALL PRBNTX t NPC5 , X X5 , Y Y5 , A NGP5 ,XHK 5 ) 

ANGP6=0. 

RFA0(NREA0,5020INPC6 
IF(NPC6 .LT. 3IG0T0 130 
REA0(NREA0,5030r (XX6( II ,I=1,NPC6I 
REAOINREA0,503OMYY6< II,I = 1,NPC6I 
UPITF(NURITE,6090) 

00 60 1=1 ,NPC6 

60 NPITE(NURITE,60 30I < X X 6 ( I 1 , YY6 d 1 1 

call PR3NTX«NPC6,XX6,YY6,AN5R6,X1X6I 

AN6R7=R5.0 

RrAO(NREA0,5070 INPC7 

IF1NPC7 .EO. 0 .AND. MSPL .tO. IIGOTO 130 

IF1NPC7 .EO. OIGO TO 80 

IFINPC7 .LT. 31G0T0 130 

READ INRE A 0,50 30 I IXX7( II,I = 1,NPC7) 

READ (NREA0,503DMYY7(II,I = 1,NPC7I 
URITEINURITE, 6100) 

00 70 I=1,NPC7 

70 UPlTE(NURITE,6030)IXX7(II,YY7ini 
UP1TEINURITF,6110)ANGP7 
CALL PR9NTX(NPC7,XX7,YY7 ,AnSR7,XNK7) 

80 CONTINUE 


A NC R 8 = 0. 

PEA0(NREA0,5020INPC8 

IFINPC8 .EO. 0 .ANO. PSBL .EO. 1150 TO 130 

IFINPCB .EQ. OIGO TO 85 

IFINPCS .LT. 3)50 TO 130 

READ INRE A 0,50 30 MXXP(I1,I = 1,NPC0I 

RFAO(NREAO,5O30l (YYB( I I ,1 = 1 ,NPC8» 

URITEINURITE ,61151 

00 8? I=1,NPC8 

82 URITE(NURITE,6D30)(XX8dl,YV8dl) 

CALL PR3MTX(NPC8,XXfl,YY8,ANGR8,X5X8l 
85 CONTINUE 



non 


«NEB9=0. 

RE«0CNBE«0»5020IN9C9 

IFINPC9 .EQ. 0 .«N0. HTC .EO. IISOTO 130 
IFINPC9 .EB. DIED TD lOD 
IFINPC9 .LT. 3IE0T0 130 
READINREAD*503aifIX9(I>,Irl,NPC9l 
REA0{NRE«D,503DMrV9III,Irl,NPC9l 
URITEINHRITE,E1?0I 

00 90 I=1,NPC9 

90 URI1E(NURITE,E030ICXX9CII,vr9(in 

GILL PRBnTXCNPC9,XX9.Tr9TlNER9*X<IK9> 

100 CONTINUE 

INBR10=0. 

REID<NRtAO,5a?0)NPC10 

IFINPCID .EO. 0 LIND. HTC .EO. 1 .AND. NCLC .EQ. IIEOTO 130 
IMNPCIO .EO. OIEO TO 120 
IFINPCIO .LT. 3IG0T0 130 
READINREI 0,5030 1 ( X XIO 1 1 1 .1 r 1 ,NPC1 0 1 
READiNREAD,50301 (YVIO (ll.Irt.NPClOl 
URITECNHRITE, 61301 
DC 110 I=1,MPC1Q 

110 URlTEtNURITE,6030l fXXlOrn ,vyiOfl>> 

CALL PRaNTX(NPC10.XX10,YY10,AMER10.XMKlDl 
120 CONTINUE 
EO TO 140 

130 URlTEtNyRITE, 61401 
GO TO 2000 
140 CONTINUE 

THE PROERAH ITERATIONS ARE CARRIED OUT TO RELITIYE ACCURACIES SPECIFIED 

BY eight CONYERGENCE factors IOENOTEO by CFXXXI. except FOR CFFLOU, 
these factors are OEFTNEO as ABSMOLD VALUE-NEU VALUEI/INEU VALUE) I. 

C CFT2 - CONVERGENCE FACTOR FOR STATIC TEHP. AT STATION 2 

C CFTS - CONVERGENCE FACTOR FOR STATIC TEHP. AT STATION 5 

C CFT5T - CONVERGENCE FACTOR FOR TOTAL TEHP. AT STATION 5 

C CEP5T - CONVERGENCE FACTOR FOR TOTAL PRESS. AT STATION 5 

C CFP4S - CONVERGENCE FACTOR FOR STATIC PRESS. BETNEEN STATIONS 4 AND 5 

C CFMCH - CONVERGENCE FACTOR FOR HACH NUHBER BETWEEN STATIONS 4 AND 5 

C CFFLOU- CONVERGENCE FACTOR FOR TOTAL INFLOW AND OUTFLOW 

C IOEFINEO as ABSMINFLOW-OUTFLOWI/TSMALLER of the two FLOWS)) 

C CFTWO - CONVERGENCE FACTOR FOR NETAL OUTER WALL TEMP. 

CFT2:.001 

CFTSr.OOI 

CFT5T:.03I 

CFPSTr.OOl 

crPAsr.ooi 

CFNCH:.OOI 

CFFL0U:.001 

CFTWO^.OOl 

C 

C SET DEFAULT VALUES 

C 

lUNTSrO 

ICTRrn 

MTC = C 

MSBLrO 

KELC-D 

RGAS-S3.35 

1 OA = 0 

DO IRE 1=1,50 
HFC4(I):i.O 
145 HFC45<I)=1.0 
150 CONTINUE 

lOArlOA * 1 

REA0(NREAD,0ATT,END=2000) 

WRITE (NWR1TE,6150)I0A 
IFdUNTS .EO. 0)G0 TO 160 
WRl TE (NWRITE.61 60) 

GO TO 170 
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160 WRITEINWIMTE, 61701 
170 CONTINUE 

IFirUNTS .EO. 0I60TO 180 
WRITE(NHRITE,6IB0IR6*S 
60 TO 190 

180 URITEINURITE,61901RG«6 
190 CONTINUE 

IFIICTR .EO. OlOHSrO.O 
IFIICTR .EO. 1160 TO 200 
60 TO 210 

200 UR1TE(NWRITE,6200I0NE 
210 CONTINUE 

IFINTC .EO. 1 .«N0. XCLC .EO. 1160 TO 220 
60 TO 2J0 

220 URITEINWRITE. 62101 
230 CONTINUE 

IFCHTC .EO. 0160 TO 290 
60 TO 2S0 

290 URITEINWRITE, 62201 
250 CONTINUE 

IFtHSBL .EO. 0160 TO 260 
60 TO 270 

260 HRITEtNURITE, 62301 
270 CONTINUE 

IFIIUNTS .EO. 1160 TO 280 
URITE(NURITE,6290INIR 
6C TO 290 

280 URlTEfNURITE,62S01NIR 
290 CONTINUE 
C 

c convert INPUT UNITS I ENGLISH OR SII TO WORKING ENGLISH UNITS 

C 

0HG:0H6*3.I 9 159/30. 

IFdUNTS .EO. OlTTrTT ♦ 960. 

IFIIUNTS .EO. 1 ITTrTT*9./5. 

IFIIUNTS .EO. 0IRS«S=R6*S*32.I79 
IFIIUNTS .EO. 1 19S«$:RG»S*5.980 

00 310 1=1, NIR 
XLOOIfII = T*UIIII/OIII I 

UR1TEINURITE,626011,MIHPRIH,DIUI,T*UIII»,*L0DIII),HSP1III,*IHPI1 
9l,RllIi,PITIII 
IFIICTR .EO. OIRIIIKO. 

IFIIUNTS .EO. 0160 TO 300 
R]III=RI 111/25. 9 
01111=01111/25.9 
TAUIIII=T«UIfIl/25.9 
HSPHII=HSPl'II 1/25.9 
XI«PIIl=XIMPIII/25.9 
PIT I IT = P1T 1 11*1.950377 
300 CONTINUE 

PITI1I=PITIII«199. 

*IHPIII=FLO*TINIH'»RII I) *3. 19 1 6* 1 01 T 1 1 /2 . 0 l**2/l 9 9 . 

01111 = 01111 / 12 . 

T«UIIII=TAUIIII/12. 

XinPIIl=XinPiIl/I2. 

HSPllII=HSPlin/l?. 

IFIICTR .EO. 1IR1III=RIIII/I2. 

310 CONTINUE 

IFIIUNTS .EO. 1160 TO 320 
WRITCINURITE,6270INFCR 
GO TO 330 

320 URITEINWRITE, 6280INFCR 
330 CONTINUE 

DO 370 I=1,NFCR 

IFIKCLC .EO. 0ITAUCII)=0. 

IFIKCLG .EO. I .AMO. TAUCfll .EO. O.OIGO TO 390 
GO TO 350 

390 URITEINWRITE. 62901 
60 TO 1000 
350 CONTINUE 

XLF C 111 = 1 TAOl II 1/SINI ALPHA IT 1/57.295781 
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XLFCPCtH = (TMHIJ»TAUC«I)>/SIMt»LPHMI)/57.2957B> 

XLFCC(I»rCT»UCJI) » /SI N ( «LPH» ( I )V5 7 . 79578 I 
XtODFC(n=XLFC(n/DFC(II 
XLODXX(nrXLFCPC(I)/OFCm 
IF(HSBL .ro. 0)R0V6(I>r0. 

IFCMSSL .EO. D»R092G( I iro. 

WRl TE(NWRITE,63D0I I ,N FCHPR til ,OFC ( I) , T*U ( I ) , TAUC ( I I .XtOOXX 1 1 I ,HSP5 
♦ til .ALPHA 1 1) .BET At II , ROVG ( I > , RO V7G 1 1 ) ,R<» tl) ,P6II ) 

IF( ICTR .E3. OIRR 1 1 l:D, 

TTtlUNTS .EO. DIGO TO 360 
RMt I l-Ril ( 1 1 //B.R 
DFCII>;3FCtI)/25.!| 

A5t I IxAStl )/|?.54l **2 
TAU(I)rTAU(Il/25.4 
TAUCtllrTAUCtll /25.4 
XlFC(I)rxLFCtII/25.4 
XLFCPCt I » ^XLFCPCtl 1/2 5.4 
XLFCC(M = XLFCCtI)/25.4 
HSP5(I)rHSP5(I)/25.4 
H GO tIl:HG 0111*0. 176228 
HPl tllrHSKIlPO. 176228 
TPSG<Il-THSGtll*9./5.-460. 

P6tl)-P6tll*l. 450377 
R0VG(I>rR0WGtI)/4.8824276 
R0V2G(I l:R0V2Gf I I /I .4881639 
360 CONTINUE 

P6t I |rP6(I 1*144. 

P5( I 1=P6I II 

A 0UT(IlrFL0ATtNFCHPRlIll*3.14lb*t0FC(Il/2.Q 1**2/144. 

A5(I)rA5(Il/144. 

DFC(ll:DFC(ll/12. 

TAUII1=TAUCI1/12. 

TAUC(IlrlAUCtIl/I2. 

XlFCtIlrXLFC(II/I2. 

XLFCPCIIlzXLFCPCtl 1/12. 

XLrCCtII=XLFCC(TI/12. 

HSP5tIl:HSP5tII/I2. 

IFtlCTR .eO> 11R4(II=R4(1I/12. 

370 CONTINUE 

IFtlCTR .EO. IIGO TO 420 

tTHf FOLLOHING CALCULATIONS ARE FOR NO CENTRIFUGAL EFFECTSI 

FIND PITHIN AND P5NAX tNINIHUN SUPPLT PRESSURE AND MAXIHUN FILH COOLING 

BACK PRESSURE 1-6ET INITIAL GUESS FOR PLENUM TOTAL PRESSURE tP3TI 

DO 380 I^I.NIR 
PlTHLOrPlTtll 
IFtI .EO. IlPlTHINrPlTHLD 
IFtPlTHLO .LI. PITHIN IPlTUNrPlTHLO 
380 CONTINUE 

OP 390 Irl.NFCR 
PEH0LD=P5tII 

IFtI .ED. I IPSNAXrPSHOLO 
IFtP5H0LD .GT. P5NAX 1 P5MAXrP5H0LD 
390 CONTINUE 

CHECK THAT PITKIN IS GREATER THAN P5NAX 

IFtPlTHIN .LE. P5NAXIG0 TO 403 
GO TO 410 

400 URITEtNWRITE.63101 
GO TO 1030 
410 CONTINUE 

P3TKXX=P1TMIN 

P3TKNNrP5MAX 

P3TrtP3TKXX ♦ P3TKNN1/2. 

60 TO 500 
420 CONTINUE 
C 

C «THE FOLLOHTNG CALCULATIONS ARE FOR CENTRIFUGAL EFFECTSI 

C find RIKN AND R4MM ILOUEST RADIUS FOR SUPPLY HOLES ANO FC HOLESl AS NELL 



C US THrtR CORRESPONOINE INOEICS fIHLO «ND JHLD I , DES I6N t T INS THE LOWEST 

C RIDIUS BT RHN. CRLCULRTE THE HIGHEST RHD LObEST SLLOUtBLE PRESSURES IN 

C THE PLENUn «T RHN WHICH PRECLUDE REVERSE FLOW IP3THXT >ND P3THNNI. GET «N 

C INITIRL PLENUH PRESSURE PROFILE MSSUNE T EOUILS TTI. 

C 

DO RIO 1=1. NIR 
RlHLO=RltII 

IFII .EO. 1>R1NN=RIHL0 
IFfI .EO. 1IIHLD=I 
IFCRIHLD .LT. RIHNIR1)1N=R1HL0 
IFIRIHLO .LT. RlHNUHLO = I 
430 CONTINUE 
RHN=R1NN 

DO 440 J=1.NFCR 
R4HL0=R4I Jl 

IFiJ .EO. 1IR4HN=R4HL0 
IFIJ .EO. 1IJHLD=J 
IF(R4HLD .LT. RRHNIR4HN=R4HL0 
IFtRRHLO .LT. R4HRIJHL0=J 
440 CONTINUE 

IFCR4HN .LT. RHNI4NN=R4HN 
P3TMXX=PlTtrHL0l 

DO 450 1=1, NIR 

PH0L0=P1 TII I*?.?!*!** »OHG*OMG*(RRN*RHN-R1 fl »*RH I J J/(?.«R5#S*TT ) I 
IFtPHOLO .LT. P3TR*XIP3THXX=PH0L0 
450 CONTINUE 

P3THNN=P6« JHLOl 

00 460 J=1,NFCR 

PHOLD=PGI Jl*2.7183**( 0n6«0HS*(RHR«RHN-R4 IU>*R4(J>)/(7.*RE»S*TT)I 
IFtPHOLO .ST. P3TRNN1 P3THNN=PH0L0 
460 CONTINUE 

IF(P3THXX .LT. P3THNNTG0 TO 490 
P?TMNR=(P3TMXX*P3THNN J/2. 

00 470 1=1, NIR 

P3TIR(II=P3TMNR*2.7183«*l0HG*0NS*tRIfI»*Rl(I)-RRN9R>1N>/{2.*RG*S*TT 
* I ) 

470 CONTINUE 


DO 480 J=1,NFCR 

PITFCRt Jl = P3THNR*?.7l 83**f ONGRORSURRt Jl *R4( J>-RHN*94N »/t?.*RG»S*T 
♦ Til 

480 CONTINUE 
GO TO 500 

490 WRITEtNWRITE, 63101 
60 TO 1000 
500 CONTINUE 
C 

C THE FLOW IS SOLVED *5 FOLLOWS - • PRESSURE OR PRESSURE DISTRIBUTION 

C (P3T OR PITIRtll I PITFCRtJI FOR NO CENTRIFUGSL *N0 CENTRIFUGXL EFFECTS 

C RESPECTIVELY! IS *5501*10 IN THE »LENUH *ND THE INFLOW XND OUTFLOW *RE 

C C*LCUL»TED FOR TH*T PRESSURE OR PRESSURE DISTRIBUTION. IHE ASSUMED 

C PRESSURE OR PRESSURE DISTRIBUTION IS THEN ADJUSTED TO E0USLI2E' THE 

C INFLOW AND OUTFLOW 

C 

C IJ IS THE COUNTER FOR THE OVERALL FLOW ITERATIONS 

C 


I J=0 

SID CONTINUE 
IJ=I J*1 


C 

C ASSUME ORIFICE TOTAL PRESSURE EQUALS SU»PLY TOTAL PRESSURE (POTtlll 

C AND THE ORIFICE STATIC PRESSURE EQUALS THE PLENUM TOTAL PRESSURE 

C IP3T OR P3TIR(II1 

C 


DO 56D 1=1, NIR 
IFllCTR .EO. DIP3TIP( I |rP3T 
P?T 1 1 ) = PI T ( II 
T2I (II=TT 
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DO 550 Ilrl,15 
P?tl )=P3TIRfI J 

IFllI ,ra. l)T2tIlr0.95rj*T2T(H 
T20-T?(TI 

CALL AIRPRPtT2D,IUNTS, 

*G .GMl .SNIDS tGPI ,GP tD2,G0G11.GNlD2,XMU,PRN.XKA,CP) 

T2HLD:T2(II 

TrRM:CI.3-:«P?«n/®2T« 1 1 l**6»1tOSI 

IFdCR'l .LT. a.O»TER?<=0.0 

V2< I>=SORT( t2.3*6*R6AS*T2r(I>/6'11l*TERH» 

XP2 < I l:V2« I )/S9RT(G*RGAS«T2T< n*CP2«I )/P2T(I M*»SH1DGI 

IFtXM?CII .SE. l.OJGO TO 520 

T?t I irI2T II »/»l .0*GH1D2*X«2C1I*X«*2( I > > 

GO TO 530 
520 XM2IH-1.D 

T2tl)d2TIIl/M.0»6HI02» 

V2( I l-SORT I6*RSAS*T2( 1 1 I 
P2(II-‘»2Tf I >/GPI0?**GD6Hl 
530 CONTINUE 

XM0=XM2II I 
NC = 5 

CALL SPLINE ( NC , NPC 5 ,X X 5 ,T T5 . XHD , AN5R5 ,X«K5,C00) 

COI (I»=COO 

RH02CT>rP2(n/(R6AS*T2(in 
T2CNVG:ABSCT2HL0-T2«1 H/T2fl» 

IFIT2CNVS ,LE. CFT2JG0 TO 560 
IFtll .EQ. 15160 TO 5<t0 
GO TO 550 

540 WRITE (NWRITE .63201 1 
550 CONTINUE 

560 RFJ2(I»-RH02(I»*V2(I» ♦0ICI»*COim2XNU 

GET AVERAGE VALUES FOR IHPIMGENENT H CALCULATIONS 

XXIHP=0.a 

XDI=a.O 

XHSPI=0.0 

XRH02rO.O 

XCOIrO.O 

XV2ro.O 

00 570 r=l,NIR 

XXIHP: XXIHP ♦ XINPIl J/FLOATINIRI 
XDI= XDI ♦ 0I(I1/FL0AT<MIRI 
XHSPl-XHSPl « HSPHIl /FLOATINIRI 
XRH02- XRH02 ♦ RH02 fl 1/FLO AT ( NIR » 

XCOI-XCOI ♦ COKII/FLOATINIRI 
570 XV2r XV2 • V2 1 1 1 /FlOA T (NiR 1 

CALCULATE INFLOW ILBN/HRl 

WIMPT-O.O 

DO 580 1=1, NIR 
C0IFC=C0ICI1 

WIHP(II=CDIFC*AIHP(II*RH02III*V2IT1*32.I74*360D. 

WIMPT=WIHPT * WINPIII 
580 CONTINUE 
C 

C 

C CALCULATE VELOCITT OUT THE FILN COOLING HOLE. ITERATE FOB P5T. 

C 

C K IS THE COUNTER FOR THE OVERALL FILH COOLING FLOW ITERATIONS 

C 

DO 760 X=I,15 
00 590 Irl.NFCR 
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JCV (I>:3 
XHBUlra. 

590 CONTINUE 

DO 670 I=1,NFCR 

IFCICTR ,E0. 0»PJIFCR«U=P3T 

IFIJCV(I) .EO. nSO TO 670 

IFIK .EO. 1 IT5TII7=TT«ITNSBtn*‘l60.-TTI*0.50 

IF(HTC .EO. 0IT5T4I»:TT 

T5T0L0(ri=T5T(U 

DO 650 11 = 1 ,15 

IFIII .EO. I IT5CII=0.95*T5TCU 
T50=T5( 1 1 

CALL AIRPRP(T5D.IUNT$, 

*6,6M1 ,SHIDG ,GP1 ,GP102,GDGMI,6H10?,*MU,PRN,XKA,CP) 

T5HLD=T5( II 

TFlPSTtU .LE. P5f 11 I P5T( I » = P5eil*l .001 
00 620 KK = 1 .15 

KCNW6 IS THE COUNTER FOR THE IHOlVtOUAL FILH COOLING ROW FLOW ITERATIONS 

KCNVGTI 1 :KK 
P5eTlrP6(Il 

IFIHK .EO. 1 IPSTd 1=P JTFCRdl 
PST0LD=P5T(T» 

TFRMz{l.Q-(P5d »/»5TI 1 I I*»6H1DGI 
IFITERN .LT. O.OITERHZO.O 

W5(IJrS0RT((2.0*G»RGAS*T5TTII/SHll*TERMl 
XMS ( IlzWSd 1/S0RT( G*RGAS*T5TT n*(P5 d l/P5Tdl)**5HnGl 
IFlXMSd) .GE. 1.01 GO TO 600 
TS(I)Tl5rdlFd.0*GHlD2«XN5dl*Xi5<in 
GO TO 613 
600 XP5dl:l.0 

T5tIlrT5T<ll/d.0»GMlD2I 
VS(llrSORT(G*RGAS*TS(Ill 
PSddPSTTI I/GP10 2**G0GH1 
610 CONTINUE 

XKOrXHSd 1 

NC = 6 

CALL SPLINE CNC,NPC6,XX6.rY6,X«T0, A W3R6,XHK6,XKT0) 
xktzxmto 

PST d i: ( P3TFCR( 1 1 ♦ P 5 ( 1 1 * XX T 1 / f 1 .0 ♦ XKTl 
PSTNEWrPST { I I 

PSTCVd IzASSlPSTNEW-PSTOLOl/PSTNEW 
IF(P5TCV(I1 .LE. CFP5TI60 TO 633 
620 CONTINUE 

HRITECNWRITF, 633011 
630 CONTINUE 

PH0Sdl=P5(Il/(RSAS*TSdll 
TSCNVGzA3S(T5HL0-T5dlt/T5dl 
IF(TSCNV3 .LE. CFT5130 TO 660 
IF( II .EO. 151G0 TO 6<I0 
GO TO 650 

ORQ WPITETNWRITE ,63101 I 
6S0 CONTINUE 
663 continue 

XXKTdlrXKT 

R0W2CnizRH05dl*V5dl*V5UI 
IFIMSPL .EO. niR0V2G( 11=1.0 

R0V2RTdl=ROV2Cdl»32.17it*36 03;»3503./R0W76(II 
IFIMSPL .EO. 31R0\/2RT d 1=3.0 
RPW2R = R3W2RT< 1 1 
XBE TA=3ETA ( I ) 

IFihSPL .EO. 0150 TO 665 


NC= 7 
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C*LL SPLINE»MC,N»C7,XX7,Yy7,R0tf??,»'IGR7,XMi<7,rD03l 
NC=e 

CALL SPLINEJNC.NPCS.XXS.YYS.XBeTA.AMGPS.XKKA.PTCGP) 

665 CONTINUE 

IFtHSBL .EO. OJCOOO-1.0 
IFJNSBL .EQ. DIRTCOP=l.n 
RTCRMirPTCOR 
CCFCII|rCDOD*RTCOR 
IF(«S8L .EO. OlROYGtl l^l.a 

ROVRAT(nrRH05(n*V5(I>*32.17«l*360a./R0VGn) 

IMNSRL .EQ. OIROYRAT ( ri:a.n 
»EJ5«I)rPH051I)*W5U)*DFCU)*CD*^C»I j/.XIU 
670 CONTINUE 

CALCULATE TOTAL AMO STATIC PRESSURE AND TEMPERATURE AT THE ENTRANCE 

OF THE FILM COOLING HOLE. 

DO 730 I:I,NFCR 
IF(JCV(I» .EO. 1»G0 TO 730 

IF(K .EO. nTAT(I)rTH(THSG«I>*RSD.-TTI*0.?D 
IFTMTC .EQ. 0)TAT«I):TT 

00 710 Il::l ,15 

IFlIl .EO. 1 IPNTd 1 XPSTI 1 1 *1 .075 
IFIII .EO. IIP<4«I)rP51I)*l.D2 
IFtll .EO. 1 1 WN ( 1 1 TVS I n*D.9B 
IFIII .EO. 1 I T4 ( I ) rTS ( H*D.?9 
P45Ir(P4Ttl).P5T(T 1 1/7.0 
P45T (P4 ( I 1 *P5(I 1 1/2.3 
W95- (V4(I 1«V5(I I 1/2.0 
T45- (T4 ( 1 l«T5(I I 1 /2.3 
P45HLn-P45 

CALL AIRPRPT I45.IUNTS. 

• 6 ,6M1 ,54103 ,GP1 ,QP102 ,G03M1 ,3M1D?,XM J,PPN,)(« A ,:p ) 

RH045RP45/(R&AS*I451 

RFYN4 5:RH045*V4 5*DFC(I)*C0FCni/<MU 

1F(REYN45 .LI. 2530.IFRFCrl6.3/REYN45 

IF(REYN45 .GE. 2 53 0 . 1 FRF C: 1 .4 22 6E -5 * RE Y N 4 F, * ] , D7 S D 9 

IF(0EYN45 .GE. 4 03 0 . 1 FRFC r 0 .395 3/ RE YN 4 5 * *0 . ?64 7 

OELPTtFRFC*XLFCPCT I l»BH045*N45**2/nrC( I) » 7.0 1 

P4T ( I irPSI ( 1 1 ♦ OELPI 

PIO^PSI ( I 1 

P INrP4T( 1 1 

A OLD -1.3 

anew:!. 3 

lonsti I 

Iiorisi ( I 1 

IIN:I4I(II 
XMOLOrXMS ( I ) 

CALL mnEa(:fmCH,PT3,pIM,A0L0,ANEW,I3,TTP,TIN,»m3L3,I'JNTS, 

♦XMNEH.INl 

IFIXMNE/ .LI. 1.3IG0 TO 690 
IFtXMMEW .GE. 1 .01 X“NEWrl .0 
T0:T4T( I 1 

call AIRPRPI TO.IUNTS, 

♦ GAH,6''1,3M106,GPI,GP107,60Gmi,34107,X4U,PRN,X4A,C'’) 

or 6flp J-l,13 

TN:T4T| II /( 1 .0 ♦ 3M1371 

CALL AIR»RP(TN,IUNTS, 

♦GTST,GM1,6M1DG,GP1,G’1D7,G0S41,54TD?.x>1U,fRN,xfA,:fi 

GCVG:A9S(GTST-3AM1/5TST 
IFIGCVS -LE. P.DOllGC 79 693 
GAMtGTST 
6R0 CONTIMJ- 

CONTI SUE 
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C»LL m^PSPdO.IUNTS, 

#G ,C-“1 ,3hIOG ,6PJ ,GP 10? ,GOGHl ,GH10?,X'«U,PON,i<'< A ,C»> 
P‘i(II=.“i4Tni/(1.0*G>11D?*X"‘iei)*Xi«III (••GOG'^l 

VU(IlrS3^U(?.Q»G*RGAt*TMT(II/G*«ll*ll.D-tP'ItI)/»ilT(I))*»GM10G)» 

RH0Am:Pq(I>/(R5AS*Ti|(I>) 

PM5N-l»i|(I)*P5(n)/?.a 
PASCNVrjjsiPASHLO-PRSMl/PRSN 
IFCPM^CNV .LE. CFPA51G0 TO 7?Q 
IF(II .E3. 15)30 TO TOO 
6P TO 71 D 

TOD WRITE tlWRI TE ,t J5Q) I 
710 CONTINJE 
7?D CONTINUE 

RAM 5 IT I rSEYF'ilS 

I''(PTC .E3. 0130 TO 770 

TC-TT-USO. 

FCHSPrnSPST 1 I 
FCHOrOF C ( I ) 

H'PrFCHSP/FCHD 
ZFC:XLOOFC( I I 
HOrHSOm 
HlrHGl t 1 I 
XlLODrxxlWP/XOI 
HFCTR-HFCM ( T 1 
HHFCTR=-IFCM5(T) 

OAU:TAU( II 
0»U?rTA JCt 1 1 
XLMrXLF; ( 1 1 
XLCrxLFcC ( 1 1 
AOSrAOUT ( I I 

IF(K .E3. 1 1 WOUT ( I 1 rw I wpt/(«-lOATT WFCRII 
WFCP:WOUT(I) 

FCBLRri WOUT ( I )/A5( I I ) 

IG;T**$GU 1 

CALL TMET0«IJ.TC,FCHSP,FCH0,H0,H1,XIL0D,xRHO2,XV2,XLM.,XLC, 
»XWSPI»HFCTR,HHFCTR ,XCOI,OAU,?FC,WFCR.A0 5 ,FCBLR,TG,HSP, I UMTS, 

*ETA,ICO,TCIN,TWI,TWO,KL’*,AKK,HO,HH,CFT«0,KCLC,AKCTtOAU?,TWIFfTCIF, 
*WPC9,WP ei3.AWGR9,aWSR tOtXX9,XXI3,ryP.yY10.XWK9,X1K10> 

TUTTIIITCIN « AbO. 

TETTIIITCO ♦ AbO. 

T*1I IIIITWI 
TPO « I i:i WO 
XFTAlTlrETA 
XHD ( I l:no 
XHHUliHH 
XXAKMHI=AKFI 
XXAKCnidAKCT 
XKLMIIlrxLJI 
TW2 ( IldWIF 
TC2«I):TCIF 
770 CONTINUE 
C 

C CALCULATE OUTFLOW <L3**/HR) 

C 

UOUTTrO. 

00 7A0 IZI.NFCR 

UOUT(T|rCOFCtl)*AOUT«H*RH05III*¥5«I)*3?.17A*3bOO. 

WOUTT^WOUTT ♦ WOUT «1) 

740 CONTINUE 
C 

C CHECK THAT T5T HAS CONVERSED 

C 

JCVT:n 

00 750 Irl.NFCR 
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n n o n 


T5TFTO:»3S(T5T(I)-T5I0L0(in/T5I«n 
ir(T5TFT3 .LE. CFFST J JCV (I »-l 
733 JCVTrJCVT ♦ JCVU» 

1F(JCVT .EO. NFCRIGO TO 7»0 
760 continue 

WRITE(NWRITEt6360)IJ 
7«0 CONTINUE 
C 

C COMPARE HEIGHT FL3US AND ADJUST P3T TO BALANCE THEM 

C 

IFIWOUTT .GT. WIM»T)FLOFCrn»OUTT-KIMPT) /UIMBT 
IFiWIMPT .oT. JOUTT »F LOFC:: I WIMPT*UOUT T ( /UOUT T 
IF(FL0FC .LE. CFFLOWIGO to 860 
IF(IJ .SE. ?5)G0 TO 8SC 
IFIICTR .EO. UGO TO 790 
C 

C (THESE CALCULATIONS APE FOR NO CENTRIFUGAL EFFECTS) 

C 

IFLUOUTI .GT. WIMOT)P3TmxK:P 3T 
IFIUIMPT .GT. N0UTT)P3TMNN:;P3T 
P3Tr(P3TMX* ♦ P3T>1NNl/2. 

GO TO 51 3 
790 CONTINUE 
C 

C (THES'' CALCULATIONS ARE FOR CENTRIFUGAL EFFECTS) 

C 

IF(WOUTT .GT. NIHPT )P3TMXX = P3TMN9 
IFLUIMPT .GT. «0UTT)P3T)«NMrP3TMNR 
P3TMN0:(P3TMXX ♦ P3TMNNJ/2. 

T8TAV-3. 

00 BOO j:l,NFCR 
TRTAWrTNTAV ♦ TNTLJl 
800 CONTINUE 

T9TAV=T8TAV/FL0AT( NFCR ) 

XTi|TAV:T»TAV-<l60. 

TITAVrUNTAV ♦ TT)/?. 

IF(MTC .£0. Q)I3TA9rTT 

ESTABLISH P3T AT THE IMPINGEMENT AND FILM COOLING BOW RAQIl ANO CHECK 

that the new PRESSURE DISTRIBUTION DOES NOT CAUSE INFLOW 

DO 810 l:l,NIR 

P3T 1R( I ) -P3TMNRP2.71B 3**L0MG*0M5*( R1 I I)*RltI)-RKN*RMN)/(?.»RGAS*TT 
*TAV » ) 

810 continue 

00 830 KNrl.lO 
JPVFLTro 
P3TFK=0. 

DO 820 J=I,NFCR 
JRVFL-O 

P3TFCR( J)=P3TMNR*?.7183**<0M6P0Mj*CRmj)*p8( J)-RKN*RMN)/(?.»R6AS*T 
93TAVI) 

IFtP3TFCR(J) .LT. P5 ( J ) ) JRWFL - 1 
JRtfFLTrjRVFLT«JRVFL 

IFIJRWFL .EO. 1 )P3THL0=P3TMNR*(1.» CP5( J) -PSTFCRC J) »/P5( J) ) 

IFLJRVFL .EO. I .AND. P3THLD .GT. P 3TFK I P3TFK rP3 T HLO 
8?0 CONTINUE 

IFIJRVFLT .EO. 0)S0 TO 8<l0 
IFtJRWFLT .GT. 0 |P3TMNR:P3TFK 
IF1P3TMMR .GT. P3TMXX )P3TNNR:P3T>)XX 
830 CONTINUE 

URITETNWRITE, 63101 
60 TO 1000 ^ 

890 continue 
GO TO SIO 

650 URITETNURITE, 63701 
GO TO 1000 
860 CONTINUE 
C 

C DATA OUTPUT DATA OUTPUT DATA OUTPUT DATA OUTPUT 

C 


60 



URlTE(NURITE,b3B0IIJ 


irtlUNTS .EO. IISO TO 870 
HRl TE f NWaiTE, 63901 UINPT 
URITECNURITE. 64001 
60 TO 880 

870 «IRPT=8I«PT*0. 45359 

WRITE (NURITE. 64 101 UIttPT 
URITEtNURITE, 64701 
880 CONTINUE 

DO 910 131, NIR 

IFIIUNTS .EO. 0160 TO 890 

PIT III3P1TC1I44.78803E-3 

P2TIII = P7Tm*4.79803E-3 

P2III3P71II44.78803E-3 

T?Ttll:T2T(I|45./9. 

T2III=T?III*5./9. 

UIHPin3UI)1PCII*a.4S3 59 
GO TO 900 
890 CONTINUE 

PlTfIl3PlTltl/144. 

P2T m = P2T(II/144. 

P2(I|3P2in/t44. 

T2TII J3T2TIII-46Q. 

72(1)372(11-460. 

900 CONTINUE 

URITE(NURI,TE,64 30)I,P1T(I), P2 ( 1 1 .XH^^( 1 1 ,T2T ( 1 1 , T 2 ( I I , tf I C I I , CDI I 

♦ I ) 

910 CONTINUE 

irilUNIS .EO. 1160 TO 920 
WRITE (NtlRITE,6440IU0UTT 
WRITE(NWRITC,645Q) 

60 TO 930 

920 WOUTT3WOUTT40.4S3S9 

WRITE(NWRITC,6460)W0UTT 
WRITE(NWRITE, 64701 
930 CONTINUE 

00 960 Irl.NFCR 

IFIMTC .EO. 0ITC2(II3TT-460. 

IFIIUNTS .EO. 0160 TO 940 
P4T(I|3P4Tm*4.7S80 3E-3 
P4( 1 13P Cl ( 11*11.788031-3 
T4T(1I3T4T(I)*5./9. 

T4(IlrT4(II*5./9. 

P5T (I)3P5T(H*4.78803E-3 
PS(1)3P5(II*4.T6803E-3 
P6( I )3Pb(I 1*4 .7B803E-3 
T5T ( I )3T5T( I )*5 .79. 

T?.( I )3T5( I )*5./9. 

TC2 (I)3(TC2( I I ♦ 460.)*5./9. 

IFIKCLC .EO. 0)TC2(I)30. 

PFCR3P3TFCR f II*4.T880 3E-3 
U0UT(TI3W0UT(II«0. 45359 
GO TO 950 
940 CONTINUE 

P4T(I )3P4T( I )/144. 

P4 1 I 13P4 ( I >/144 . 

T4T(I)3T4T(I)-460. 

T4( I >374 ( 1 I -460. 

P5T (I )3P5T ( 1 1 7144. 

P5( I >3P5 I I > 7144. 

P6(I )3P6( 1 17144. 

T5T I I I3T5T( n -460. 

T5(I)3T5(l)-463. 

IFIKCLC -EO. 0)TC2(II30. 

PFCR3P3TFCRI1)7144 . 

950 CONTINUE 

WPITEINW. 4 HE, 6490)1, PFCR,P4(I), *14111, T4t(ll, 1411), P5TII), “5(11, 

• XH5(1I,T5T(I1,T5(II,TC2(1I,W0UTII),**KTH1,C0FC(I),9TCIM1),R0V4AT( 

*1),R0V2RT(I),KCNV5(I) 



n r> n 


960 COMTIMUr 

IFIHTC .EO. 0»&0 TO inoo 
WR1TE(NWRITE, 64951 
IFCIUVTS .EO. 1)53 TO 970 
WPITE»I9M91TE,6490) 

GO TO 930 
970 CONTINUE 

«RITE(NW»ITE, 65001 
980 CONTINUE 

DO 995 In.NFCB 

IFdUNTS .EO. 1)50 TO 985 

»5)I)r65U)9144. 

GO TO 990 
985 continue 

H GO (I)=HSO(I >*5.67446 

HEl II)=HGim*5. 67446 

XHH(I>:KHH(I)*5.67446 

XHDtI)'KHD(I)*5. 67446 

*5II)rA5(I)*9?9.0304 

TNS6I I) - CTNSSII ) <460. > *5. /9. 

TNO(I)-(TMO«I1 ♦ 460. 1*5. /9. 

TW2m:fTH2<I) • 460.)*5./9. 

IFIKCLC .EO. 0)TW?)I)-D. 

T)«I fl)dTNI(I) » 460. 195. /9. 

XXAKHtll rXXAKNd )*0.0I7296 
XX«KCTd)=XX*KCT(I l«0.0t7?96 

9»0 MRITE (NU9ITE,65I0) I.HEOd) .HGId) ,XHH tl ) ,XHD ( I ) , HF C4 5d I •HFC4(1 1 , 
«ASd>.THSG(I>,IMOd),TU?(I>.THICI>,XXAKN(I>,XXA<CTlI),XET«d),KKLH 
* d> 

995 CONTINUE 
1000 CONTINUE 

RETURN VARIABLES TO ORIGINAL TN^UT UNITS 

ONErOMGtIQ./I. 14159 
IFdUNTS .EO. OlTTrTT-460. 

IFdUNTS .£0. llTT:TI*5./9. 

IFdUNTS .EO. 0IRSAS = RGAS/3?.I74 
IFdUNTS .EO. 1)R5AS:R6A5/5.990 

DO 1020 1 = 1, NIH 
IFdUNTS .EO. 1)50 TO 1010 
OTd)=OI d)*12. 

TAUT (I)=TAUId)*12. 

H5P1 (I)=HSPld)*12. 

X1MP<I)=XINP(I)*12. 

Rid )=R1 d)*12. 

GO TO 1020 
1010 CONTINUE 

RId>=RIdl*304.8 
DId>=DId)*304.8 
TAUI (I)=TAUId)*304 .3 
HSPl{I)=HSPld)*304.8 
XlHPTI)=XINPd)*304.8 
1020 CONTINUE 

DC 1040 I=I,NFCR 

IFdUNTS .EO. 1150 TO 1030 

0FC«I)=0FCI1)*12. 

HSP5(II=HSP5(I»*12. 

TAUCI)=TAUf 1)*I2. 

TAUCIII=TAUCIII*12. 

R4d)=R4d)«12. 

IFtHTC .EO. DTASTl )=ASdl*144. 

GO TO 1040 
1030 CONTINUE 

R4d>=R4d)«304.8 
OrCII) = OFCd)*304.B 
TAUd) = TAUdl*304.8 
TAUCIII=TAUC{II*304.8 
HSP5d)=HSP5(I»*304.8 



ROVEI n=R0VEIII*<>.e82H276 
ROV2EU l=ROVZEI I 1*1 .«881b39 
IFIMTC .EO. D»»5fT»:*SIIM9?9.030<» 

10«t0 CONTINUE 
C 

C format STATEMENTS 

C 

5010 FORMATClbASl 

5020 F0RMAT(I2I 

5030 FORMATISFIO.OI 

6010 FORMAT! IHI ,// ,lbA5 ,// I 

6020 FORMAT! /, IX ,57H 

• ,/,5X ,93HIN!>UT POINTS FOR COOEANI GAMMA VERSUS T ARE.F.lZX.l 

•HX,9X,lHr ,/ 1 

6030 FORMATtSX .2F10.ll I 

6090 FORMAT!/, IX ,57H 

• ,/,5X,it7HIN»UT POINTS FOR COOLANT VISCOSITY VERSUS T ARE./.l 

72X, 1HX.9X ,1HV ./ > 

6050 F0RMATI5X,FI0.<l,2X,E12.<ll 

6060 FORMAT!/, IX ,57H 

• ,/,5X,51HIN»UT POINTS FOR COOLANT SPECIFIC HEAT VERSUS T ARE 

* ,/, I2X, IHX ,9X ,1HY,/ ! 

6070 FORMAT! /, IX ,57H 

• ,/,5X ,58HIN?UT POINTS FOR COOLANT THERMAL CONDUCTIVITY VERSU 

*S T ARE,/,12X.1HX,9X,IHY,/1 

6OR0 FORMAT!/, IX, 57H 

• ,/,5X,A9HINPUT POINTS FOR IMP. OISCH. COEFF. VERSUS M2 ARE,/ 

*,12X,1HX,9X,1HY,/I 

6090 FORMAT! /, IX, 57H 

« ,/,5X ,67HlNOUT POINTS FOR FILM COOLING TOT. PRESS. LOSS COEF 

♦F. VERSUS M5 AR E , / , 1 2 X , 1 HX , 9X , 1 HY , / > 

6in0 FORMAT! /, 1 X ,57H 

« ,/,5X,ii9HIN“UT points FOR FILM COOLING VERSUS R0V2» ARE,/ 

*,12X,1HX,9X,1HY,/I 

6113 FORMAT!/, 5X,16MR0TATI0N ANGLE : , F 1 0 . 3 ,2 X , 7HD E GRE E S , / ) 

6115 FORMAT!/, IX, 57H 

, ,/,5X,3SHINOUT POINTS FOR RTCOR VERSUS BETA ARE , /, 1 2X , IHX ,9X 

»,1HY,/1 

6120 FORMAT! /, 1 X ,S7H 

• ,/,SX,MRHIN»UT POINTS FOR MEIAL CONDUCTIVITY VERSUS T ARE./, 

»1 ?X , 1HX,9X , 1 HY,/ 1 

6133 FORMAT!/, IX, 57H 

• ,/,SX,50HlN°UT POINTS FOR COATING CONDUCTIVITT VERSUS T APE, 

•/,12X,1HX,9X,1HT,/1 

61<iD FORMAT!/, SX.SRHCASE ABOOTEO - A REOUIRED CURVE VAS NOT INPUT OR WA 
•S SPECIFIED BY LESS IHAN 3 POINTS./) 

6150 F0RMAT!1H1,10X,28H OUTPUT FOR C HA MB ER, 1 5 , 1 3H ,/l 

6160 FORMAT!/, 5X,18HSI SYSTEM OF UNITS! 

6173 FORMAT!/, 5X,23HENGLISH SYSTEM OF UNITS) 

61B0 FORMAT! /,5X .2IHC03LANT GAS CONST A NT r , n , F 1 D . 3 ,2X , 8 H J/ ! KG -< ) ) 

6193 FORMAT!/, 5X,21HC00LANT GAS CONS T A NTr , IX , F 1 0 . 3 ,2X , 1 6H I F T-L B F ) / ! L BM- 
*R 1 1 

52'’3 FORMAT( /,5x,63HTHIS CASE INCLUDES CENTRIFUGAL EFFECTS. ROTATIONAL 
*SPEED E3UALS,F10.?,2X ,NHPPM.| 

6213 FORMAT! /.SX.BBHTHIS CASE INCLUDES A THERMAL BARRIER COATING) 

6223 FORMAT! /, 5X ,7BHTHIS CASE TS FLOV ANALYSIS ONLY ANO INCLUDES NO mft 
*AL TEMPERATURE CALCULATIONS) 

6233 FORMAT!/, 5X,36HFHrS CASE HAS NO MAIN STREAM BLOVING) 

62M3 FORMAT!///, 1 X , I 5 , 2 X , 2 5 H ROW S OF IMPINGEMENT HOLES,//, 5X, THROW, 2X, 
♦5HH0LFS,2X, 1 3H0I AMEIER ! IN) ,9X .RHVALL ,8X ,3HL/0,9X ,9HH0LE,5X,nHIMP 
»I NGEMENT ,6X ,2HR1 ,9X,3HPlT,/,33X,?HTHICKNE5S,16X,7HSPACINo,RX,8HPIS 
*TA.NC£,6X,itH(lNI,6X,6H!PSIA),//) 

6253 FORMAT!///, 1 X , I 5 , 2X , 2 SHROUS OF IMPINGEMENT HOLES,//, 5X,3HROU,?X. 
»SHH0LES,2X,13H3IAirTEP IMMl,<tX,RHVAlL,BX,3HL/0,9X,ilHH0LE,5X,nHlMP 
*IN'GEMENI,6X,2HR1,9X,3HP1T,/.33.X,9HTHICKNE5.S,16X,7 HSPACING,iiX,SHPIS 
*TANCE,6X,MH!MMI,5X,9HIN/CM»*2),//) 

6263 FORMAT! 5X,I3,i»X,t3,RX,F7.9,6X,F7.3,5X,F7.3,5X,F7.3.5X,F7.3,5X,F7. 
•3,5X,FB.3) 

6273 FORMAT!///, 1 X , I 5 , 2 X , 2bHR 0 VS OF fjlm COOLING HOLES,//, 1X,3HR0V,2T 
* ,5HH0LES,?X ,13H0IAMETER !INI,7X,9HIHICKNE5S,7X,3HL/0,9X,RHH0LE,SX, 
*5HALPma , 

» 5X,iH3ETA,7X,5HRM3VG,7V,6HRH0V?S,9X.2Hi5t|,6X,2HO6,/,33X,I5HWALL- 



* C0«TING,2X,TH<T0TlkLl .6X, 

«7HSPACING,3X,5Hf DrGi.SX,5H(0ESI ,2X,1 <iH(LBH/FT** 2*HRI ,]X,14H(LBM/FT 
**HR**2I .IX.RHf IN>,3X,6H«PS1»1 ,//l 

6280 F0RHAT(//2, 1 X, IS , 2X, 26HR0US OF FILM COOLING HOLES, 2/, 1X,3HR0U,2X 

♦ ,5HH0LES,2X,13H01AHFIER I PHI ,7X ,9HTH1 CKNESS ,7X ,3HL /O , 9X , RHHOLE . 5 X , 
ASHALPHA , 

* ,SX,4H3ETA,7X,5-IRHOV6,7X,6HRHOV2S,9X,?HP4.6X,2HP6,/,30X,1SHRALL- 

♦ C0ATING.2X .7HIT0TAL ) .6X, 

*7HSPACING,3X, SHIOEGI, 5X,SHI0EGI,2X,12HIK G/MAA2AHRI ,2X,l?HrX3/M*HR* 

• «2) ,<tX,4H(MMI ,lX,9H(N/CH*A2l,//> 

6290 FORHAK /,5X,6?HCASE ABORTEO - COATING RAS SPECIFIED BUT NOT COATIN 
*G THICXNESSI 

6300 FORMAT! 1 X ,I 3 ,4 X ,I 3 ,4 X , F 7.4 ,SX,F7 .3 , 3X , F 7 . 3 , 3X , F 7 . 3,5 X , F 6 . 3 , 4 X , F 6 . 

*3,3X,F6.3,2X,E12.5,F12.5,4X,F7.3,?X,F8.31 
6310 FORMAT!//, 5X,66HCASE ABORTEO - THE SPECIFIED PRESSURES KILL RESULT 

* IN REVERSE FLOW, //I 

6320 FORMAT!/, 5X,67HWARNING - T2 HAS NOT CONVERGED IN 15 ITERATIONS FOR 

♦ IMPINGEMENT ROW, 151 

6330 FORMAT!/, 3X,59HUARNING-P5T HAS NOT CONVERGED IN 15 ITERATIONS FOR 
*F.C. ROW, 15) 

6340 FORMAT!/, 5X,68HWARNIN6 - T5 HAS NOT CONVERGED IN 15 ITERATIONS FOR 
A film cooling row, 15) 

6350 FORMAT!/, 5X,1IIHHARNING - THE AVERAGE PRESSURE BETWEEN STATIONS 4 
AAND 5 HAS NOT CONVERGED IN 15 ITERATIONS FOR FILM COOLING ROW ,15) 
6360 FORMAT!/, 5X,73HWARNING - T5T HAS NOT CONVERGED IN l5 ITERATIONS IN 
AOVERALL FLOW I TER A T ION , I 5 I 

6370 FORMAT!/, 5X,70HIMPINGEMENT AND FILM COOLING FLOWS HAVE NOT CONVERG 
AEO IN 25 ITERATIONS) 

6380 FORMAT! //,5X,52HIMPINGEHENT AND FILM COOLING FLOWS HAVE CONVERGED 
AIN, I 5,2X ,1 SHOWER ALL I TERATl ONS , // // ) 

6390 FORMAT! TOX , I3HINFL0W EQUALS, F9 . 3, 2X , 6HL B M/HR ,/// ) 

6400 F0RMAT!3X,3HIMP,2X,6HPSPLYT,6X, 2HP2 , 5X , 2HM2 ,6X , 3 H T 2 T , 5 X , 

A2HT2,7X,4HWIMP,3X,SHC0IMP,/,3X,3HR0W,2X,6H(PS1A ),21X, 

A 3H!F) ,1 2X,8H!L3M/HR) ,/T 

6410 FORMAT! lOX, 13HINFL0W EOUALS.FR. 3, 2X , SHK G /HR , / // ) 

6420 F0RMAT(3X,3HIMP,2X,6HPSPLTT,6X, 2HP2 , 5X , 2 HM? , 6 X , 3 H T 2 T , 5 X , 

a2HT2,TX,4HW1MP,3X,5HC01)1P,/,3X,3HR0W,1X,9HIN/CMAA2) ,19X, 

A 3H!K),13X,7H!KG/HR),/) 

6430 F0RMAT12X,I3,2X,FT.3,3X, F7.3,1X,F5.3,4X,F5.0,2X,E5.'3,3X,FT 

A.3,2X,F5.3) 

6440 FORMAT!///, IOX,14HOUTFLOW E OUA L S , F 9 . 3 , 2 X ,6HL 3 M/ HR , // / > 

6 4 50 FORMAT! I X , 2HFC , 4 X , 3HP 3 T , 7X , 2HP4 , S X , 2HM4 ,3X,3HT4T,4X,2HT4,1X,1H/,3X 
A , 3HP5T,bX ,2HP5,5X,2HM5,3X,3HT5T,4X,2HT5, 1X,1 H/, 5HTCTIF,4X,4HW0U 
AT,3X,2HKT,4X,?HRT,5X,2HRT,4X, 

A 4HRHOV,4X,5HRHOVSO,IX,4HITRS,/,1X,3HROW,?X, 

A6H!PSIA),17X,3H!F),TX,1H/,1X,6H(PSIA),18X,3H(F),7X,IH/,1X,3H!F),3X 
a,8H!LRM/HR),13X,4HC0RR,3X,5HRATI0,3X,5HRATI0,/) 

646D FORMAT!///, 10X,14HOUTFlOW E OUA LS, F 9 . 3 , 2 X ,5 HM G /H R , // / ) 

6470 FORMAT! I X , 2HFC , 4 X , 3HP 3 T ,TX , 2HP4 , 5 X , 2HM4 , 3X , 3H T4 T , 4 X , 2H I 4 , 1 X , 1 H/ , 3X 
a,3HP5T,6X,2HP5,6X,2HM5,3X,3HT5T,4X,2HT5,1X,IH/, 5HTCTIF,4X,4HW0U 
aT,3X,7HKT,4X,2HRT,5X,2HRT,4X, 

A 4HRHOW,4X,5HRHOWSQ,1 X ,4HI I RS,/, I X ,3HR0W, 1 X , 

a9H(N/CMAA2),15X,3H!K),7X,1H/, 9H!N/CMAA7),16X,3HIK), 7X,1H/,IX,3 

AH(K),3X,7H(KG/HRI,14X,4HCORR,3X,5HRATIO,3X,5HRATIO,/) 

6480 F0RMAT!IX,I2,1X,F4.3, IX,F8.3,IX,F5.3,1X,F5.D,1X,F5.D, 1H/, F«. 

A3,lX,rB.3,lX,F5,3,lX,r5.3,lx,F5.3, IH/, F5.0 , 1 X , F T . 3 , 1 X , F 6 . 3 , 
A1X,F6.3,1X,F6.3,IX,F7.3,1X,FT.3,3X,I2) 

64 9D FCRMATI///,2X,2HFC,2X , 26HHE A T -TRANSF E R -C OFF F IC I E N T 5 , 2X , 1 3HH -MOD -F A 
aCTORS,4X,6HCOOLEO,4X,3HGAS,8X,I6HWALL-IEMPERATURE,7X,17HAVG.-THFRM 
A.-C0ND.,5X,3HET A ,7X,3HITR, 

A /, 1X,3HR0W ,2X,3HHGD,3X,3HH51 ,2X ,7HF C -H OL E , 1 X , 4 H I M? 5 , 5 X , 7H 

AFC-H0LE,2X,4HIMPG,5X,4HAREA, 5X,4HTEMP,4X,7H0LITSIDE,2X,THTNTFACE,2 
ax,6HINSI0E,3X,5HMETAL,3X, THCO A T TNG , 3X , 9 H ! TCO-TC) /, 

A /,13X,16H!BTU/FTaa2*HPaF),22x,7H!INA*2),4X,3HIF),6X, 

A3H!F),5X,3H!F),5X,3H!F), 7X , 1 3H I RTU/F TAH RAF ),4X,8M!TW0-TC),/) 

6495 FORMAT!//, 1CX,21hHEAT TRANSFER PESULIS) 

6500 FORMAT! ///,2X,?HFC,2X ,26HHEAT-TRANSPER-C0EFFICIENT5,2X,1 3HH-M0D-FA 
ACTORS, 4X,6HCOOLEO,4X,3HG AS, ?X,)5HWALL-TEMPERATURE,7X, I 7HAW3. -THERM 
a.-COND.,5X,3HEIA,7X,3HITR, 

A /, IX ,3HR0W ,2X ,3HMGD,3X,3HHG1 ,2X , 7HF C -HDL E , 1 X , 4 H I M PG , 5 X , 7H 

aFC-H0LE.2X,4HIMPG.5X,4HAREA, 5X,4HTEM‘>,4X,7H0UTSI0E,2X,THTNTFACE,2 



o o non 


*X,6H INSIDE, IX, 5HMET«L,3X,7hC0«TTN6,3X,9H(TC0-TC)/, 

« /,lJX,I6H(J/(»«*«i*SEC*Kn,22X,7H(CH**2),9X,3H(K),6X, 

«3H(K ) ,SX, 3H(K ) ,5X,3H«K ) ,6X, 14HIJ/CC>**SEC*K )) , <1 X , 8H ( T W 0- TC ) / I 
ftSIO rORMATI lX,I2,IX,'^S.0,lX,Eb.D,IX,P5.0,lX,Ffc.0,?X,F5.3,3X,F5.3,3X,F7 
».3,5X,FS.0,4X,F5.7,IX,FS.3.3x,F5.0,<IX,F6.7,<lX,F5.3,3X,F7.i|,bX,12 1 
GP TO 153 
2033 S’OP 
END 


SUSPOUTINE TNEI3 


SUfcROUTINE TMET0IIJ,TC,FCHSP,FCHD,H0,H1,XIL0D,XRH02,XV2,XLM,XLC, 

• XHSPl,HFCTR,mFCTP,XC0I,QAU,7FC,WFCR,A05,FC8LR,TG,HSP. lUNT 

♦ S,ETA,TC3,TriN,TWI,Ti(0,KLH,»KN,HD,HM,CFT«0,XCLC,»KCT,0»iJ2,TWlF,TCI 
•E, 

*NFC9,NPC10,ANGR9,AN5R10,XX9,XX13,yT9,TY10,XNK9,XHX10) 

DIMENSION CMAT{?9,2‘:i ,CN(?4> 

DIMENSION XHHU2‘»»,XMK?(2»I,XNK3(2iU,XMK9(2‘»),XMi<9l29),XN<iai29> 
DIMENSION XXl(24l,XX2(24».XXIt24»,XX4f29),XX9(29),XXlDI2‘l) 
DIMENSION YYI(24l,YY2(24I,YYSe24),»Y4r24),YY9(24l,YY13l24I 
COMMON NPCI,NPC?,NPC3,NPC4,ANGR1,ANGR2,ANGR3,ANSR4,XX1,XX?,XX3,XX4 
*,YYl,YY2,YY3,YY4,XHKI,XM<?,XMK3,XM<4,NREAn,NWRIIE 

TPIF-T5-IC 

TCINrTC*T0IF*.2D 

IWI:IC»TDIF«.25 

IClF:TC«rDIF*.3D 

IwrF:TC»TQIP*.35 

TC0rTC*T3IF*.43 

IW0:IC«T3IF*.4S 

OP 70 KLMrl ,15 
TOLOnwO 

IFILMrO.S*(TUI«tCI»4tiO. 


call AIRPRPdFlLM.IUNTS, 

♦G.GMl ,&M10G,GPl ,G»1D2 ,60GM 1 ,6M 192 ,XMU ,PR N ,XK A ,CP I 
-arrival velocity '■ACTOR 


IFIXILOD .LT. 3.780IFACVA:I.0 

IFIXILOD .GT. 3.T80IFACVA:-.D39I93*XILOD»XILOO 
*♦ .93715 

IFIXILOD .GT. 7.590IFACVAr0.902597*XIL0D*XIL0D 
** 1.46D5I9 

IFIXILOD .GT. 14.34»FACVAr0.901390*XIL0D*XIL0D 
•« 1.IP5854 

RENArrACYA»XRH02*XV2*XHSPI/XMU 

PENA:.RENAPXC0I 


.051495*XIL00 

.107042AXIL0D 

.0SD825*XIL0D 


-GARDON AND COBONPJE IMPINGEMENT CORRELATION 

HPrlD.28t,*XKA/XHS»lT«RENA**D.b25 

H0:H0*HFCTR 

AREARrl.3-3.1RI59/I4.C«HSP*HS?» 

HC^AREARPHD 

TOArO.S*! TC0«TCIN) *460. 

IFIKClC .EO. II TCA:D.5*(TCIF»TCINI*46D. 

CALL AIRPRPITCA.I'JNTS, 

*G ,GM1 ,GMI 0G,GP1 .GP1D2 ,GDGM I .GM102 , XMU.PRN ,XK A ,CP I 

TR:lITilO«T«lJ/2.0)»itbD. 

IFIKCLC .EO. 1 >TR=I ITMIF^TDI »22.0»«4b0. 

TWAVrTR 



irduVTS .CO. 1 >T JAV: T9*5./9. 

NC-9 

tall ^PLINEINC,NPC9,XX9,YY9,TWAV,AXS99,XMK9,A!<1) 

miUNTS .EO. 1 ) A<MrAK«*57.817& 

I'^CKCIC .E3. 1 1 TCTAV= t TUir*TwO» /?. » Ofcri. 

IFCKCLC .EO. 1 .AMD. lUNTS .ED. I ) T C T A V : T CT A v *5 . / 9 . 
I*’(KCLC .EO. 0)AKCT = 3. 

IF(KCLC .EO. 0>G0 TO 10 
1C 


CALL 'iOLIMF(NC,MPClD,YX10,YY13,TCTAV,ANGR10,X**.<lD,Ai<CT) 

IFlIUMTS .EO. 1 )AXCT:AKCT*57.R17b 
10 CONTIMUE 

U=HC*OAU/AKM 

TP; ( TCA/T!J)**Q. 1 3 

H IP THE FILM COOLIMG HOLE CALCULATED FROM T.°. DAYEY C D 7 7£ L A T I ''M 

CC‘EF:(3.035*XKA*T7/FCH0I*«1.0/T'’C1**D.? 

PFH-(WFCR*FCHril/ ( A05*XMU*J2.17T9«7b0n. ) 
HH-C0FF*«REH**3.8»«(PRN**0.M) 

HM=HH*T.imS9*FCH0*FCH0*ZFC/(FCHSP»FCHSP*nAUl 

HM:HM«HHFCTP 

OA:HM*OAU*OAU/AKM 

TFtKCLC .EO. nXLTOTiXLM*XLC 

IFtKCLC .EO. nTRr(((TCO«TCTF)/?.*3SD.I/irTUO«TOIF)/2.*4b3. ))♦*.)' 
IF(KCLC .EO. IIGO TO ?0 
GO TO 30 

?0 TCCAV:(TCIF«TCO)/?. 

call AIRPRP(TCCAV,IUMTS, 

*G,6M1 ,GM1 DS.GPl ,6P1D2 ,G0GM1 . 3M1 0 2 , X MU, P P N , XK A , C » ) 

REh: (WFCRAFCHO l / (A05*XMU*32.17 39*3600. ) 

30 CONTINUE 

IFIKCLC .EQ. llHH2rO.OM3*XKA*<REH**D.8>*IPRN**0.it)*TR*tF:HD**0.?l* 
* <XLTOT**0.8-XLM**0.8)/(FCHO*(XLTOT-XLM) ) 

IFIKCLC .EO. 1 lHM2 = HH2*3.mi59*FCHD*XLC/ (FCHSP*F:HSP*0AU7) 

IFIKCLC .EO. I ) HM?:HH2*HHFCTR 
IFIKCLC .EQ. llDA?rHM2*0AU2*0AUZ7AKCT 
TCAP:n.5*(TC0«TCl»M6a. 

CALL AIRPRP«TCAO,IUNTS, 

*G,GMI,GM10G.GP1,GP1D2 ,G0GM1 ,GMI02,XMU,PPN,XKA,CP> 
BFTA=DAU*HM/(FC8LR*CPI 

IFIKCLC .EO. 1 )BrTA230AU2*HM/(FC3LR*CP) 
ROOTr(BETA*BFTA*<I.O*DAI**0.5 

IFIKCLC .EO. llROOT2::(BETA2*BETA2«it.O*OA2l**0.5 
Alr-0.5*IBETA»RD0T) 

A2::-0.5*IBETA-R00T1 

IFIKCLC .EO. 1)AL1=-0.5*IBETA?»R90T21 
IFIKCLC .EO. 1 »AL2:-O.5*IBETA2-R00T2) 

IFIKCLC .EO. IIGO TO 40 
DFNr|U-A?)*EXPIAl)-tU-An*EXPIA21 
C2^ IU-A2)/DFN 
C3rlAI-Ul/0EN 

ETArC2*tl.0-Al*A170Al*EXPIAl»»C3*ll.D-A?*A2/DA)*EXPIA2) 

IFIETA .GE. l.OlETArO.9999 
DELHG-HO-Hl 

TWO = TG-IITG-TCI*IE TA*FCBLR*CP»0.0-ETA1*OELHGIMIHD-ETA*DFLH6»ETA* 
*FCBLR*CP> 

TNEWrTWO 

TCO-ETA* miO-TC 1«TC 

TCIN=IC2*I 1 .D-A1*A 1 /DA l«CJ*I 1 .0-A2*A2/DA 1) • I TUO-TC I *T C 
THIrlC?*C31*IT«0-TC)*TC 
GO TO 50 
40 CONTINUE 

CMATI1,1I=U-A1 



CHJmi.ZIrU-*? 

CM*TIi,31:0. 

CHAT(],«|rO. 

CN«T(l,5|rO. 

CH*T(i,6»=0. 

Cn«TC?,l |=( AKM/AKCT>«(0«U?/0*Ut*kl*eXP(A 1) 
CNAT(?,;i=(AKn/AKCTI*(0AU2/0AUt*A7*EXP(A2> 

CHAT(?, JI=0. 

CHATC2,A|r-ALl 

CMATI?,5|r-AL? 

CMATI2,5I=D. 

CPATn,H=EXP<Al> 

CMATl3,?J=tXP(A?) 

CMATC3,3I=-I. 

CMAT«3,<»lr-l. 

CHAT<3,5)r-l. 

CPAT»3,blrO. 

CPAT(a, 1 >rC1.0-AI«Al/DAI*EXPIAll 
CPAT<4.2|r(I.Q-A?*A?/0AI*EXP(A2l 
CKAT(4,3>=-1. 

CHAT Cu.m r-« 1 .0 -AL 1*A L 1 /OA 2 I 
CHAT(4,5»r-»1.0-AL2*AL2/0A2» 

CHAT<t|,b»=Q. 

CMAT(5,U=0. 

CHAT(5,2)=0. 

C HA T ( S , 3 ) = 1 . 

CHATt5,M |rEXP(ALl » 

CHAT(5,5lrEXP«AL2) 

CHAT(S,bl=l. 

CALL XHrxSLfS.CHAT.CNI 

C2:CN(1 I 
C?:CNJ21 
C<t = CN»3| 

CS=CNf4) 

CbrCN(S| 

ETA=Ca»C 5*< l.a-ALl*ALl/PA2»*EXO(«Ln«Cb»n.D-AL2«»L2/0A2)*E*P( AL2» 

IF<ETA .GE. 1.0>ETAr0.999P 

0ELHG=H3-H1 

TWOrTG-((TG-TCI»(rTA*ECBL9*CP*tI.O-ETAI*DFLHGI)/(HO-ETA*3ELHG»ETA« 

*FCBLP*CP1 

TNEU=TXO 

TC0:ETA*«T«0-TCMrC 

TVIF5(C2*EXP(AlI»C3*EXP»A2)T*flU0-TC»»TC 

TriFr(C2*U.O-Al««l/DAt*EXPrAIM:i*tl.O-A?*A2/DA)«EXPtA2n*(TMO-TC 
* 1«TC 

TCIN:{C2*(1.0-Al*Al/0AMC3*U.a-A2*A2/0«n*tTyD-T!:)*TC 

TWlr(c2«C3)*(TM0-TCI*TC 

■=0 continue 

TU0CVGrABS(TNE«-T3L0IFTNEW 
TF(KLH ,E0. IlGO TO 70 

iFtTuncxG .le. cftnogo to PO 

IF(KLH .EO. I5IG0 TO 60 
GO TO 70 

6L WRITE<NURITE,6u0nj 
7t continue 
8U continue 

IFIKCLC .EO. OITWIFrO. 

IFCKCLC .EO. 0)TCIF:0. 

brc F0RHAU/,5X,93HWARNING - OUTER W«LL IcHPEPATU = E h'S NOT C f N V t pG f 
♦IN IS ITERATIONS IN OVERALL FLOW J TER A T ION , I S t 
RE turn 
END 


SUB ROUT I NE HNE V 


Sl'BROUT INE XNE V< C^MCH ,PTO.PTN.AOLO,«NEV.TO,T1C',tIN,X'<CL'-',iL'VlS, 
♦XHNEU.TN) 



u u u 


OinCNSlON XX] (24 1 ,XX2(2M,XX3f2«l*XX4l2<» I 
OIHENSION YYl(24l,Yr2(24l,VV3«24l.yr4(24 > 

OIHENSION XnKl(24l ,XHK2(24I ,XNK3(24I,XRK4(24 I 

COHHON NPCl ,NPC2tNPC3,NPC<.IIN6ftl,«NGR2,«NBR3,XNSR4,XXI tXX?,XX3,XX4 
•,YY1,YY2,YY3,YY4,XHKI ,XRK2 ,X«R3,XBX 4 ,MRE RD.NMRITE 

CALL A1RPRP(T0,1UNTS, 

• 60,EHI,SNI06,6PI ,SPI02.E06NI,SRI02,rr(U,PRN,XKA,CPI 

FOR THE first ITERATION EVALUATE BARHA AT THE GIVEN TOTAL TEHP. 

AND LET THE FIRST GUESS AT XHNEU BE XHOLO 
C 

CALL AIRPRPITTN.IONTS, 

*6N,GHI,GHID6t6PI tSPI02,G0GHt ,GN102,XHU,PRN,XKA,CPI 

XHM=XHOLO 
DO 10 1 = 1 ,25 

PATG=(PT0/PTN»*»A0L0/ANEVI*SQRTIITN/TT0>*S0RT(S0/GN» 

XHFCN=1.D * ( (GN-1 .01 /2.0I*XNN*XHN 
XMFCO-l.O ♦ ((GO-1 .01/2. OI«XHOLD»XHOLO 
P0WNr(GN*1.0l/(2.0*(GN-l.DII 
POUO=(60*l.rM/(2.O*(G0-l.Oll 

XNUH-XHN-PATG*XM0L0«( XHFCN>*4P0MN/t I XMFCO)**POVO I 

DNM=l.a-PAT6*(GN*l .01 •XHOIO*XHM*XHFCN**( (-GN«3.ai /(2.0*(6N-1 .0111/ 
* (2.0*XHFC04«P0V0) 

XHHLO-XMN 

XHN=XHN-*NUH/ONM 

TN=TTN/(1 .0*( (GN-1 .0) /2.0I*XHN*X4NI 
CALL AIRPRPdN.IUVTS, 

*GN,GMl,SM10G,GPt.GPI02,GOGHI,6HID2.XnU,PRN,XKA,CPI 

CNVCR^ABS t XMHL0-X4N1/XMN 
IF(CNVCR .LE. CFHCHIGO TO 20 
10 CONTINUE 

MPITE (NURITE,600I 
?0 CONTINUE ' 

XNNEwrXVN 

buc FORHAT (/,EX,46HHARNING - K HAS NOT CONVERGED IN 2E ITERATIONS) 
RETURN 
END 


subroutine airpbp 


SUttPOUTINE AIRPRPI TO, lUNTS, 

*t ,GH1 ,o«10G ,GP1 ,G’102 ,G0GH1 ,GH102,X>«U,PRN,XH A,C® ) 

DIMENSION XXl(2m,YYl(24 1,XMKl(24) 

DIMENSION XX2(24),YY2< 24 ),XMK 2(241 
DIMENSION XX3(2M),YY3(24»,XHK3(24) 

DIMENSION XX4(2m,YY4(2M),XMK4(24l 

COMMON NPC1,NPC2,NPC3,NPC4,4NGR1,ANGR2,ANBR3,ANSR«,XX1,XX2,XX3,XX4 
*,YYl.Yy2,YY3,YY4,XMt(l,XMK2,XHX3,XMK4,NREAD,NVRITE 

IFdUNTS .EO. ) )TD:T3*S./9. 

NC=l 

CALL SPL1NE(NC,NPC1,XX1,YTI,T0,ANGRI,XMK1,G) 

NC-2 

CALL SPLINElNC,NP:2,XX2,Yr2,T0,ANGR?,XMX2,XMUr 
NC= 3 

CALL SPLINE (NC,NPr3,XX3,YY3,T3,4NBR3,XMK3, CP) 

NCL4 

CALL SPLINE(NC,NP:4,XY4,YY4,TD,A'IBR4,X'1K4,X<A) 

irilL'VrS .EO. 1 )XMU-XMU*. 067197 
IFlIUNIS .EO. 1 )CP=C?*0. 23901 
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o n r> 


iniUNTS .EO. I ) X<A:xh»*57.8176 
IFdUNTS .EO. 1 ) TO^^TD^'J./S. 

PPNrC'»*XMU»'»bOO./XK« 

6 P 1 rG -1 .0 
GPinG^GHl /G 
&P1 -G*l .0 
GP1D2-GO1/2.0 
GDGPjn .0/G>*1D3 
GPlD?-&«l/2.0 

THE FOLLOWING NU WftS DIHENSION OF SLUG/ T F T«SE C.I 

XHU-XMU/32.1739 

PFTURN 

END 


SUBROUTINE “RBIiTX 


SUBROUTINE PRBMTX(NPl,XK,F,ftNGR.SOLI 
C 

C THIS SUBROUTINE GENEROTES THE PR08LEH H«TRIX (H»TII,jn F«ON THE 

C INPUTEO X AND Y VALUES AND CALLS XHTSOL TO SOLVE IT 

C 

DIMENSION XK(7<t I ,F (2V I ,XKR (?AI ,FIM2(| I .SOL(?A > 

OTHENSION XXl(24UYTt(24l,XHKl(24t 
DIMENSION XX2I24 »,YY2 (24 »,XMK7C2«I 
DIMENSION XX3(24>,YY3l24l,XHK3t24l 
DIMENSION XX4(24),YY4(24I,XMK4(24) 

REAL L(24I,LR(24»,MAT(24,25) 

INTEGER OPT 

COMMON NPCI ,NPC2fNPC3 ,NPC4 , A NGR I » A NGR 2t A NGR 3 , AN 6R 4 , X X I , X X ? , X X 3 , XX U 
*,YYI,YY2,YY3fYY4,XMKl , XMK 2 , XMK 3 .X NX 4 ,NR t AD ,NWR I T E 

NrNPl-1 

NP2tn«2 

OPin 

IFIANGR .EO. D.DIOPTiD 
IFiOPT .EO. DIGO TO 20 
ANGR0TrANGR*3.14IS93/I80. 

CANrCOSI ANGPOTI 
SANrsiNCANGROTI 

00 10 Id.NPI 

XKR(II=XX(II*CAN ♦ F(I1*SAN 
FR« I »=F( IMCAN - XK(II*SAN 
ID CONTINUE 
20 CONTINUE 

DO 30 1=1, N 

IFCOPT .EO. DIGO TO 30 
LP(II=XKR{I*1I-XXR(I) 

30 L(Il = XK(I*ll-XX(n 
C 

C GET SLOPES AT THE END POINTS 

C 

YPFSTrlF(2J-F(l»»/<XX(2l-XXfIII 
YPLST = (F(NP1I-F(NI )/IXK(NPU-XN(MII 

IFIOPT .EO. 1 JTPFSTR=(FR(2»-FRIIII2(XKR«2I-XXR(11 I 
IFIOPT .EO. llYPLSTR=IFRINPII-FRINII/tXKR(NPll-XKR(NI) 

C 

c INITIALIZE THE ENTIRE MATRIX TO ZERO 

C 

00 40 1=1, NPl 
00 40 J=1,NP2 
40 MATII,JI=0. 

MAT(l,ll=L(ll/3. 
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DO r> o o 


MAT (I,?)=L(I)/b. 

MATIl,NP2»r»F(2»-F(I|)/L(n-YPFST 

PAT(NPl,N)rL(N>/6. 

MATINPI.NPI )=L(N»/3. 

MATtNPl,NP2l:YPl.ST-(F(NPl)-F(M)12L(N) 

IFtOPT .FQ. IIMAH 1 ,1 1 -LSU ) /?. 

IF(OPT ,E0. 1)MAT( 1,2 l-LRI 1 >/6. 

IFJOPT .EO. l)M«Tei,NP2)-(FP(21-FR(l ) »/LR«II-YPFSTR 
IF(OPT .EO. IIMSTCNPI ,N)rLR(N)/6. 

IFtOPT .EO. 1 )M«T( NPl ,NPU:LR( N J/3. 

IF(CPT .EO. 1 )M«T( NPl ,NP21 =YPLSTR-C FR (NPl )-FR(N) ) /LR« R) 


DO 50 Ir2,N 

IH2- 1-2 
IP1=1»1 

MAT ( I,IH1 ) IMH /6. 

MAT « 1 ,I » = (L nni ) *L t U > /3. 

MAT« I ,1P1 I=t t U fb. 

MAT( 1 ,NP2 l = «F (I PI I -F( I n/L ( I 1 - T F ( I ) -F ( I M 1 1 ) /L ( I M 1) 

IFIOPT .EO. 1 >MAT( I ,IMl JrLRUMl I /5. 

IFtOPT .EO. I IHATf I ,I »::(LR(IMI l*LRt I I >/3 . 

IFtOPT .EO. l»MATtI,IPI>:LRtII/6. 

SO IFIOPT .EO. 1 )MATtI,NP2JrtFRIIPII-FRtI) »/LRII )-tFRtTI-FRIIMl )I/LR( 
*IM1 I 


CALL XMTYSUNPl ,MA T.SOL > 


RFTURN 

EK'D 


SUBROUTINE SPLINE 


SUBROUTINE SPLI N E t NC , NPl ,XK ,F, X , A NSP , XHK N ,Y I 
C 

C THIS SUBROUTINE GIVES 4 CURVE. FIT VALUE OF Y FOR A SPECIFIED X 

XMKNI24I IS THE SOLUTION VECTOR OBTAINEO FROM THE INPUTEO X AND Y 

VALUES IN SUBROUTINE PR8MTX 

REAL MKM1,MK,LK 

DIMENSION F(2RI ,XK (2>t I ,XMKN(2<) I 
DIMENSION XKRT2A I ,FRt 2<t ) 

DIMENSION XXlf2RI,YTlt?A>,XHKt(2i|| 

DIMENSION XX2f2i|l,YY2(2>II.XMK2(2iU 
DIMENSION XX3I2R > , YY3 t 2R I , XMX 3 1 2 A I 
DIMENSION XX4(2<ll,YY4f2A>.XMK«t2AI 

COMMON NPC1,NPC2,NPC3,NPC4,ANSH1,ANGR2,AN6H3,ANSR<* ,X X 1 , XX 2 , XX 3 , XXN 
♦,YY1,TY2,TY3,YY4,XMK1,XMK2,XMX3,XMKR.NREAD,NWRITE 
INTEGER OPT 

N rN P 1 -1 
NMlrN-1 
OPTrl 

IFIAN6R .EO. O.OIOPT=D 
IFIOPT .EO. OIGO TO 15 
ANGR0T=AN6R*3.14l593/tSD. 

CAN=C0SIAN6R0TI 
SAN=SINIANGROT) 

00 ID Irl.NPl 

XKRCl)rxKtII*CAN • FIIMSAN 
FRtI>=Flll*CAN - XKII»*SAN 
10 CONTINUE 
15 CONTINUE 

FOR A GIVEN X, FIND THE XX THAT BRACKET IT AND CALCULATE GENERATED F 
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IND=0 

IF(X .CQ. XK(1IIIND=-1 

ircx .lt. xk(iiiino='2 

IMINO .LT. OiyrFm 

if(ino .lt. oigo to 80 

IFIX .ED. XK(NPIMIN0-1 
1F(X .GT. XK(NPI >IIN0:2 
IFIIND .GT. OITrFINPlI 
IFIINO .GT. OIGO TO 80 


00 30 1=2, NPI 
IND=0 

IFIX .E9. XKdDYrFdl 
IF(X .EO. XKdlieO TO dO 

IFIXKII-II .lt. X .«N0. XKdl .GT. X I GO TO 20 
GO TO 30 
20 CONTINUE 

1 PI = 1-1 


PKNlrXHXNdPl 1 
PXrXHKNd I 
XXM = XKd-l I 
XXtiXKII I 


FX = Fd» 
FKHl=Fd-l » 

LK = XKIII-XX ll-I I 


IFIOPT . 
IFIOPT . 
IFCOPT . 
IFCOPT . 
IFCOPT . 
GO TO <10 
30 CONTINUE 
<10 CONTINUE 
IFCOPT . 


EO. 

EO. 

EO. 

EO. 

EO. 


EO. 


1 1 XXH = XXR CI-1 I 
nxX = XKPC II 
1 IFXrFRCl I 
1 IFKMl =FR (1-1 1 
1 ILX=XKR( I 1-XKRCI-l 1 


OIGO TO TO 


TIL=CNKHl*CXX-XXHI**3l/Cb.*lKI*(FXNl/LX-LK*MKPl/G.)*CXX-XXNI 


V2L:(I.0/T*N(*nGR 0TII *XXM-X/SIMC»NGROTI 
IFCtU .GT. V2L1I'(0C=I 
IFCY2L .GT. YILIINOC:-! 

INOCP=INOC 

INDCPI=-INOCP 

OELXR=(XX-XXM1/10. 

XR=XXM 


00 50 1=1,30 
XPrXR»DELXR 

TrRNir.(MXMl*(XX-XRI*«31/Cfc.*LK> 
TERN2=CNK*(XR-XXN>*»3I/C6.YLXI 
TrRMI=CFKXLX-l1X*LX/6.l*(XR-XXNI 
TFRP<t = (Fxm /LK-LKFMKIC 1/6. IFCXX-XRI 
YI=TERNI*TERM2*TERM3»TERMR 
Y2= ( l.Q/TXN CANGROT n*XR-X/SIN(ftN5R0H 
IFIYI .ST. Y2)IN0C=I 
IFCY2 .GT. Y1I1N0C=-1 
CP1T=ABS(TI-Y2I/ABSCY1I 
IFCCRIT .LE. 0.0002IG0 TO 60 
IFCINOC .EO. IN0CPII0ELXR=-DELXR/1D. 
IFlINOt .EO. INOCPI)INDCP=INOCPI 
IFCINOC .EO. INOCPI )INOCPT=-INOCP 
50 continue 
60 continue 

ANGINV=-ANGROT 
SANI=SINIAN61NVI 
CANI=C0S(ANGINVI 
Y=Y1*CANI - XR*SANI 
GO TO 110 
70 CONTINUE 

TERNdCMKKMCXX-XlAYll/CG.ALKI 

TERH2=(MK«(X-XXMI*«3)/C6.»LXl 

TERP3=(FK/LK-MX*LX/6.)*CX-XXN1 

TERHA=IFKH1/LK-LX»MXH1/6.1*CXX-XI 

Y=TERNI»TERN2»TER“3*TERNA 



non no 


LE. 1160 TO no 


60 TO 110 
80 CONTINUE 


IFI 

IND 

.BE. 

-1 .AND. 

IND 

IF! 

I NO 

.EO. 

-21 GO TO 

90 

IFI 

IND 

.EO. 

2IG0 TO 

100 

GO 

TO 

no 




90 UR1TE(NURITE,600M,NC 
GO TO 110 

100 URITEINVRITE,610l)ttNC 
no CONTINUE 

600 F0RH»TI/,5X,31HM*RNING - A SPECIFIED X-VALUE (,Ft0.3,3SHI IS BELOW 

♦ THE RANGE OF INPUT TABLE, 131 

610 F0RHAT(/,SX,3IHUARNIN6 - A SPECIFIED X-VALUE (,F10.3.35HI IS ABOVE 

* THE RANGE OF INPUT TABLE, 131 
RETURN 

END 


SUBROUTINE XNTXSL 


SUBROUTINE XMTXSLINR.XNAT.SOLI 
C 

C THIS SUBROUTINE TAKES THE PR08LEK MATRIX AND SOLVES IT BT THE GAUSS- 

C JORDAN ELIMINATION METHOD 

C 

C NR IS THE number OF ROWS IN THE MATRIX lORDER OF MATRIX) 

C XMATfI,JI IS THE PROBLEM MATRIX TO BE SOLVED lINCLUOINS THE FORCING F) 

C XMAT{I,JI IS READ IN CONTTNUOUSLT BY ROWS CINCLUDING THE '"ORCING FUNCTION) 

C 

C M*T(I,J) IS THE OVERALL MATRIX OBTAINED BY ADDING THE IDENTITY MATRIX 

C TO THE PROBLEM MATRIX 

C 

C sold) IS THE SOLUTION VECTOR 

C 

DIMENSION SOLCZRi ,FCT (?<i l,XMATf 2<t,2S> 

DIMENSION XXKZB ),TYI(2<||,XMKIC?4) 

DIMENSION XX?(^C||,YV^(^4|,XMK?f?>l) 

DIMENSION XX3(2BI,YY3(?NI,XMK3I?>II 
DIMENSION XX4(24),VY4<2<l),XMKiH2it) 

COMMON N^Cl ,NPC?,NPC3 ,NPC4 , ANGRl , AN GR? , A NGR3 , ANS R4 ,XX1,XX?,XX3,XX4 
•,YY1,YY?,YY3,YY4,XMK1,XMK2.XMK3,XMK4,NREAD,NWRITE 
REAL MAT)34,49I 

NM=NR-l 

NC:NR«1 

NN:NC«1 

NLST3NC«NR 

DO 10 J=1 ,NC 
DO 10 1=1 ,NR 
ID MAT(1,J)=XMAT(I,JI 
C 

ADD THE IDENTITY MATRIX TO GET OVERALL MATRIX 

DO 30 J=NN,NLST 
DO 30 1=1, NR 
MAT ( I,JI =0. 

IF<«J-I) ,E0. «NR»niGO TO ?0 
GO TO 30 
?0 MAT(1,J)=1. 

30 CONTINUE 

make the pivot element the LARGEST ELEMENT 

NSW=0 

DO SO J=1,NR 

IF(J .ED. NRIGO TO 6D 
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nan a a n non r» o n 


or 50 lrj,NH 

IP:I«1 

IF( ABS( H«T( IP.J) I .LT. «BS(MAT( J, Jl I 150 TO SO 

NSUrNSK«I 

00 40 JS=1,NLST 

ST0R = W»T « J, JS» 

HATCJ.JSIrHXTdP.JS) 

MAT «1P,JSI = ST0R 
40 CONTINUC 
50 CONTIMUE 
60 CONTINUE 

REDUCE ELEMENTS IN PIVOT COLUMN TO 7ER0, EXCEPT PIVOT 

DO 80 J=I ,NP 
DO 70 IRri.NR 

70 FCT llR>rM«TllR,JI/MAT ( J.Jt 
FCTIJKO. 

00 80 1?ER=1,NR 
DO 80 J?ERrj,NLST 

MAT ( IZER, JZERI = MAT ( 1 2 ER , JZE P l-F CT 1 1 ZE R J * MA U J , JZE R ) 

80 CONTINUE 

GET THE OETERMINANT 

OETrl.D 
DO 90 K^IiNR 
90 0ET=DET*HAT (K ,K » 

OET=OET*< (-1.1**NS«I 

TRAP SINGULARITY 

1 SNGLrO 

IF(ABS(MAT(NR,NRII .LT. l.E-7 .AND. ABSIDETt .LT. l.E-7150 TO 100 
GO TO 110 
100 CONTINUE 
ISNSLrl 

URITE(NURITE,600I 
no CONTINUE 

DIVIDE EACH ROM BY IT’S PIVOT TO GET SOLUTION VECTOR ANO INVERSE MATRIX 

DO 120 IPIV:|,NR 

OIVrMATCIPIV.IPIVI 

00 120 JPIVri.NLST 

MATflPlV, JPIVI=MAT(IPIV,JPIVI/OI» 

IZO CONTINUE 

DO 130 10 = 1 ,NR 
130 S0L(10I=MATII0,NCI 

600 FORMAT!/, 10X,36HSIN6ULAR MATRIX IN SUBROUTINE XMTXSLI 
RETURN 
ENO 
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TABLE I. - INPUT DATA FORM 


Repeated 
for each 
input 
table 


Repeated 
for each 
chamber 


HUE 

P90JCC1 «Wt8t9 

MAtrSI 

lHt£l OF- 

SralEMCKl 

FORTRAN STATEMENT 

CC7>t8<AiaN 

1 2 3 fc « 

1 7 8 9 c 0 a 


■9 20 2> 22 23 2a|? 3 26 21 » 29 vd B 32 33 36 33 SAjlT 38 39 AO At a2|a3 AA AS A6 A7 A8 |a9 M W 32 VS 6A 

33 36 37 38 34 6oIm 62 63 6A 63 66l67 68 64 70 D 72 

73 7A 73 76 77 78 79 80 

■ 

■ 


■ 

1 




1 








1 


■ 




■ 

1 








$ 





B 





i 


r^E 













1 

1 


MMI 

IBI 












Bi 



Bil 

Bi 








bb 

■ 

BBBB 

QHAMI 1 
_$BAT 
NIR = 

ER TKPU 

T VARI 

^LES 

IN NAM 

EL 1ST 

FORM 

( 

SELECT 

ED VAR 

BH 

SHOWN 



.IUNT5 

^ 0 , .. . 



. NIHP 

R=9 .10 

.3*12. 







^B 




mp 

7 .. NFC 

ilPR= 7* 

15 . DF 

snn 

1 

0. 5* 

liTVIi 








■BBH 








■■1 






Imbi 



■p 



HHH 



imiii 



bbsBI 


■BiH 

bib 


BB 

1 

mm 

mmm 






bbbbI 


isam 

34E12 S 

bbb 

BB 

1 



BBli 



BBB 


bbbb 

■■1 

iBn 


BBB 

BB 

1 

B 

Q2EE 

BBbI 


B8 

BBB 









1 


HHH 

BBB 



bbb 

mgm 





■uni 











■■1 

!BB 













Bil 






Si 

■ 1 








3 U 13 16 17 18 

19 20 ^ 22 23 2A 

23 26 21 28 24 » 

3> 32 33 SA 33 36 

37 38 34 AO U 42 

A3 AAA3A6A7 Al 

14 30 3i 32 33 3A 

33 36 37 38 34 Aojoi 62 63 6A 63 66 

67 68 64 70 71 72 

7S lA 73 76 77 78 79 80 
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TABLE n. - TABLE INPUTS FOR FCFC PROGRAM 


Table 

Table variable, 

Correlating 


y 

parameter, 

X 

1 

Coolant specific -heat 

Coolant temperature. 


ratio, 

T„, K (°R) 

2 

Coolant viscosity, 




g/cm ■ sec (Ibm/ft- sec) 



3 

Coolant specific -heat at 




constant pressure, C 
J/g-K (BtuAbm °R) 



4 

Coolant thermal conductivity. 




kj,, J/cm-sec-K (Btu/ft- hr- °R) 



5 

Impingement -hole discharge 

Impingement -hole 


coefficient, (CD)j 

Mach number, M 2 

6 

Film -cooling -hole .total -pressure 

Film -cooling -hole 


loss coefficient, (KT) 

* ' *nmg 

Film -cooling •hole flow reducticm 

Mach number, M- 

7 

(PV2) /(pv2). 


due to main -stream -gas flow at 
S = 0°, RT 



8 

RT correction factor. 

Compoimd angle, 0, 


(RT)y(RT)^^o 

deg 


9 

Metal thermal conductivity, kj^^. 

Metal temperature. 


J/cm-sec-K (Btu/ft- hr- 

T„, K (°R) 

10 

Ceramic coating thermal conduc- 

Coating temperature. 


tlvlty, k ., J/cm-sec-K (Btu/ 
ft-hr-°R) 

T^8 f 

ct* 

(°R) 
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table m. - CHAMBER INPUT VARIABLES 


(a) Variables associated with types 
of calculations desired 


Variable 

Description 

Type^ 

lUNTS 

Input units - 0 for U. S. cxistomary 

T 



units; 1 for SI units (default = 0) 



ICTR 

Centrifugal effects - 0 to exclude; 




1 to include (default = 0) 



MTC 

Metal temperature calculations - 0 




to exclude (flow analysis only); 
1 to include (default = 0) 



KCLC 

Coating - 0 for no coating; 1 for 




coating (default = 0) 



MSBL 

Main -stream blowing - 1 for blowing; 




0 for no blowing (default = 0) 



OMG 

Blade rotative speed (default = 0. ), rpm 

R 


(b) Impingement -hole-row variables 


NIR 

NIHPR 

R1 

DI 

TAUI 

HSPl 

XIMP 

PIT 

TT 

RGAS 


Number of impingement -hole rows (=^5) I 
Number of impingement holes per row I(NIR) 
Radial location of each impingement row R(N1R) 
from shaft centerline, mm; in. 

(Input only if ICTR=1) 

Hole diameter of each impingement row, 
mm; in. 

Impingement -insert thickness at each 
row, mm; in. 

Impingement -hole spacing at each row, 
mm; in. 


Impingement distance between insert and 
shell inner surface at each row, mm; 
in. 

Supply total pressure at each impinge- 

9 ' f 

ment row, N/cm ; psia 
Coolant siipply total temperature, R 

K; °F 

Coolant gas constant (default = 53. 35), R 
JAg-K; ft-lbf/lbm-°R 


^Where 1 denotes integer; R denotes real; NIR denotes 
number of impingement rows; and NFCR denotes 
number of film -cooling rows. 





TABLE m. - Concliided. 


(c) Film -cooling -hole -row variables 


Variable 

Description 

Type^ 

NFCR 

NFCHPR 

R4 

DFC 

A5 

TAU 

TAUC 

HSP5 

HFC4 

HFC45 

ALPHA 

BETA 

HGO 

HGl 

TMSG 

P6 

ROVG 

ROV2G 

Number of film -cooling rows (^50) 

Number of film -cooling holes per row 
Radial location of each film -cooling row 
from shaft centerline, mm; in. 

(Input only if ICTR=1) 

Hole diameter of each film -cooling row, 
mm; in. 

Shell outer -surface area associated with 

2 2 

each film -cooling row, cm ; in. 

Shell metal thickness at each film -cooling 
row, mm; in. 

Coating thickness at each film -cooling row, 
mm; in. (Input only if KCLC=1) 

Film -cooling -hole spacing, mm; in. 

Local back -side impingement -heat -transfer 
correction factor. (Default=l. 0. ) 

Film -cooling -hole heat -transfer correction 
factor. (Default=l. 0. ) 

Film -cooling -hole inclination angle at each 
row (fig. 3), deg 

Film -cooling -hole compound angle at each 
row (fig. 3), deg 

Main-stream heat -transfer coefficient at 

coolant outlet temperature equal to main- 

2 

stream -gas temperature, J/(m -sec-K); 
Btu/(ft^- hr- °R) 

Main-stream heat -transfer coefficient at 
coolant outlet temperature equal to shell 
outer -surface temperature, J/(m -sec-K); 
Btu/(ft^-hr-°R) 

Main -stream -gas temperature at each film- 
cooling row, K; °F 

Main-stream -gas static pressure at each 

O 

film -cooling row, N/cm ; psia 
Main-stream -gas density times velocity, 
kg/(m^- hr); lbm/(ft^- hr). (Inpirt only if 
MSBL=1) 

Main -stream -gas density times velocity 
squared, kg/(m-hr^); lbm/(ft- hr^). 

(Input only if MSBL=1) 

I 

I 

F 

(NFCR) 

t(NFCR) 


^Where I denotes integer; R denotes real; NIR denotes 
number of impingement rows; and NFCR denotes 
number of film -cooling rows. 




TABLE rV. - INPUT LISTING FOR EXAMPLE PROBLEM 


Card 

colvunn 


Tabular 

inputs 


r 


Chamber J 
inputs I 


Title card 


>'c 


Hr VS- T- 


C_ „ vs. T„ 
p,c c 


k„ vs. T„ 
c c 


(CD)j vs. 


M 


(KT)nmg^"- “ 


RT vs. 
(RT) 


(pV^)e 




^^^/3=0° 
k„ vs. T 


(pv2)g 

vs. 


Tct 


Exampie 1 


Exampie 2 



2 EXAMPLES FOR FCFC PROGRAM 


{ 

{ 

{ 


{ 


10 


300* 

2100. 

500* 

2500* 

700. 

900* 

1100* 

1300* 

1500. 

1800* 

1.40 
1 *270 

1*386 

1*238 

1.365 

1.345 

1.329 

1.316 

1 *304 

1 *288 

300. 

500* 

700. 

1000. 

1500* 

1900. 

2500. 


1 .800E-4 

7 

2.6SQE-4 

3.350E-4 

4.200E-4 

5.400E-4 

6.300E-4 

7.600E-4 


300* 

500* 

700. 

1000. 

1500* 

1900. 

2500. 


1 *004 

7 

1.025 

1*067 

1.138 

1 *234 

1*305 

1.546 


300* 

500* 

700* 

1000. 

1500* 

1900* 

2500* 


2.510E-4 

3.849E>4 

5.062E-4 

6.862E-4 

9.414E-4 

1.172E-3 

1.736E-3 


10 

0.0 

.95 

.05 

1*0 

*20 

.30 

• 40 

.55 

.70 

*85 

*80 

.9225 

*8025 

.9225 

• 8175 

• 840 

*875 

.8975 

.91 

.92 

10 

0*0 

*95 

*05 

1.0 

*20 

*30 

.40 

.55 

.70 

*85 

*85 

*50 

.8475 

*4665 

*84 

.8275 

.805 

-750 

• 665 

.5675 

10 

0*0 
1 *0 

.01 

3.2 

.03 

*06 

.10 

.20 

.40 

.60 

0.0 
. .945 
3 

0* 

1*0 
’ o 

*20 

1*0 

45* 

2.0 

*55 

90. 

1.9 

• 6B 

*76 

.86 

.91 

.93 

T 

700. 

811. 

922. 

1033. 

1144. 

12S6. 

1367. 

18??. 

1700* 

• 2525 

• 2802 

.3111 

• 3125 

.376? 

.8116 

.8862 

• 4635 

. .576 


r 3 

\ 1033. 1811* 2387* 

^ .0131 .Omv . .0163 

^ SOflTT 

IUNTS=lf ICTR=0, MTC=1» KCLC=lt MS6L=lt RGASr287.05, 

NIR^St NIHPRrio«ls, DI: 10«0 *3048 » T AUI :1 0«0 * 6 35 , 

HSPiriO*3*81, XIMprlO*! *27, PlTr 10*R0*» . , TT-811., 

NFCR:*», NFCHPR=25*15, OFC=3*0*279q , 0*25R0t 21*0.* *5 = 25*0.9677*1* TAU=25*1 .27, 
< TAUC=2590.127, 

MSP5 = 25*2.5<» , HFCq=25*l., HFC95=25*1.* ALPHA = 40.* 38** 35*. 33.* 21*0** 
BtTA=25«0** P6=373.4* 370.8, 368.6* 364*7, 21*0.* TMS6=25*2550* * 

H60=5277., 5816*, 6384*, 6951., 21*0** H&lr3972.* 4256** 4483., 4767., 21*0., 
R0VG = 4 .364C;6* 4 * 78 1E6 , .5 • 1 07E6 * 5.59E6, 21*0., 

V. P0V2G = 3.204E12, 3.872E12* 4.439E12* 5.362E12* 21*0*S 
r lOATT 

ICTR=1, hTC=0, KCLC=0, 0M6=1682S., 

NIR=15, NIHPRz15*2, 

R1Z217.2* 219*7, 222*3, 224*8* 227.3* 229*9, 232*4, 235*0, 237*5* 240*0* 242*6 

245.1, 247.7. 250*2, 252*7, 

01 = 15*0.4318, TAUI=15«0*381, HSP 1 =1 5*3 * 8 1 * X I MP= 1 5*0* 76 2 , TT = 8ll** 

P1T=284*3, 286*4, 288*5, 290*7, 293.0, 295*2, 297*6, 299*9, 302*3, 304.8, 
307*3, 309*8, 312*4, 315*1, 317*8, 

NFCR=15, NFCHPR=15*2, 

"< R4=217*2, 219*7, 222*3, 224.8, 227*3, 229*9, 232.4, 235*0, 237*5* 240.0, 242*6 

245.1, 247*7, 250*2, 252*7, 

0FC=15*0.4572, TAu=1S*1 *016, HSP5 = 15*2 * 54 , AL PHA = 1 5*30 . , BETa = 15*0., 

P6=264*S, 265*2* 266*0, 266*8* 267*6* 268*4* 269*2, 269*9, 270.7* 271*5, 272*3 

273.1, 273*8, 274*6, 275*4, 

R0VG = 5 .089E6, 5.108E6, 5.127E6, 5.145E6, 5.163F6, 5.182E6, 5.200E6* 5.219E6, 
5.237E6, 5.256E6, 5.275E6, 5.293E6, 5.311E6, S.330E6, 5.348E6* 
R0vZG= 5.873E12* 5*902E12, 5.930E12, 5.9S8E12* 5.987E12, 6.015E12, 6*044El2, 
6.C72E12, 6.101E12, 6.129E12, 6.158E12, 6.186E12* 6.215E12, 6.243E12, 

^ 6.272E12S 
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TABLE V. - TITLE CARD AND TABULAR DATA 


OUTPUT FOR EXAMPLE PROBLEMS 


2 EXAMPLES FOR FCFC PROGRAM 


INPUT POINTS 

FOR COOLANT GAMMA VERSUS T ARE 


X 

Y 


300.0]]0 

1.4000 


soo.oooa 

1*3860 


733.0000 

1.36S0 


900*0000 

1.3>iS0 


1100.0000 

1.3290 


1333.3330 

1.3163 


1500.0000 

1*3040 


1800*0000 

1.2880 


2133.0000 

1.27Q0 


2500*0000 

1.2360 



INPUT POINTS 

FOR COOLANT VISCOSITY VERSUS T ARE 

X 

Y 

333 .3333 

. 1800-03 

500.0000 

.2650-03 

700.0000 

.3350-03 

1333 .3000 

.■1200-03 

1500.0000 

.5400-03 

1900.0300 

.6300-03 

2533.3000 

.7600-03 


INPUT POINTS 

FOR COOLANT SPECIFIC HEAT VERSUS T ARE 

X 

y 

300.0000 

1.0040 

500.0000 

1.0250 

733.0333 

1.0673 

1000.0000 

1.1363 

1500.0000 

1.2340 

1933.0000 

1. 3050 

2500.0000 

1.5480 


INPUT POINTS 

FOR COOLANT THERMAL CONDUCTIVITY VERSUS 

X 

Y 

333*0000 

.2510-03 

500*0000 

.3849-03 

700*0000 

.5062-03 

1033.0333 

.6862-03 

1500.0000 

.9414-03 

1900*0300 

.1172-02 

2533.0000 

.1736-02 


input points 

rOR IMP. OISCH. COEFF* VERSUS M2 ARE 

X 

Y 

.0330 

.8000 

• 0500 

.8025 

.2000 

.6175 

*3000 

.8400 

.4000 

.8750 

.5500 

.6975 

*7000 

*9100 

*8500 

*9200 

*9500 

.9225 

1 *0000 

.9225 


ARE 
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TABLE V. - Concluded. 


INPUT POINTS 

FOR FILM COOLING TOT. PRESS. LOSS COEFF. 

X 

¥ 

.0000 

.8500 

.0500 

.8475 

.2330 

.8403 

.3000 

.8275 

.4000 

.8050 

.5500 

.7500 

.7000 

.6650 

.8500 

.5675 

.9500 

.5000 

1.0000 

.4 665 


INPUT POINTS 

FOR FILM COOLING RT VERSUS R0V2R ARE 

K 

¥ 

.3300 

.0000 

.0100 

.2000 

.0300 

.5500 

.0600 

.6 800 

.1000 

.7600 

.2000 

.8600 

.4000 

.9100 

.6000 

.9 300 

1.0000 

.9450 

3.2000 

1.0000 

ROTATION ANGLE : 45.000 DEGREES 


INPUT POINTS 

FOR RTCOR 

VERSUS BETA ARE 


X 

¥ 



.3333 

1.0033 



45.0000 

1.0000 



90.0000 

1.0000 



INPUT POINTS 

FOR metal 

CONDUCTIVITY VERSUS T 

ARE 

X 

¥ 



700.0000 

.2525 



811.3333 

.2802 



922.0000 

.3113 



1033.0000 

.3425 



1144.0000 

.3762 



1256.0000 

.4 116 



1367.0300 

.4462 



1422.0000 

.4635 



1 700.0000 

.5780 



INPUT POINTS 

FOR COATI 

N6 CONDUCTIVITY VERSUS 

T ARE 

X 

¥ 



1333.3333 

.0131 



1811.0000 

.0149 



2367.0000 

.0163 




VERSUS MS ARE 



TABLE VI. - EXAMPLE 1 (VANE) CHAMBER OUTPUT 


OUTPUT FO» CHAHBCR 1 

SI ststep of units 

COOLANT GAS CONSTANT: 257.050 J/(K6-K> 

THIS CASe INCLUDES A IKEPNAL BARRIER COATING 


3 ROWS OF InPlNGENENt HOLES 


ROW 

HOLES 

diaheter (NMI 

WALL 

THICKNESS 

L/0 

HOLE 

SPACING 

ihpingehent 

distance 

R1 

(RNl 

PIT 

IN/CN«*2I 


15 

.3090 

.635 

2.083 

3*810 

1.270 

.000 

ADA •UDU 

2 

15 

• 30A8 

.635 

2.083 

3.810 

1.270 

• 000 

ADA .OOU 

3 

15 

.3098 

• 635 

2.083 

3.810 

1.27U 

.000 

AOA .000 


A ROWS OF FILM COOLING HOLES 


ROW 

HOLES 

OIANETER INNI 

THICKNESS 
WALL COATING 

L/0 

( totali 

HOLE 
SPAC IN6 

ALPHA 

<OES> 

BETA 

IDE61 

RKOVG 

fK6/N**2*HPI 

RN0V2G 
IK6/HAHRAA2 1 

RA 

(NH 1 

P5 

(N/CN«*2» 

, 

15 

• 279A 

1.27U 

.127 

7.779 

2.5AU 

AO. OOU 

• QUO 

.A36A0«07 

.320AU*13 

.uou 

375. AUU 

2 

15 

• 279A 

1.270 

• 127 

8.121 

2.5AO 

38.000 

• 300 

.A78 10*07 

.38720*13 

.000 

370.500 

3 

IS 

• 27WA 

1.270 

• 127 

5.717 

2.5AO 

35.000 

• QUO 

•51070*07 

.AA390* 1 3 

.000 

368. 5U0 

A 

15 

• 25A0 

1.270 

.127 

10.098 

2. SAG 

33.000 

.000 

.55900*07 

.536?0«13 

.OOU 

36A. 70U 


IHPlNGLMENf AND FILR COOLING FLOWS HAVE CONVERGED IN 8 OVERALL ITERATIONS 


INFLOW eOUALS 20.263 A6/HR 


IHP 

HOW 

PSPLTT 

TN/CNAP27 

P2 

N? 

T2T 

(K> 

T2 

WIHP 

(KG/HRI 

COIHP 

t 

AOA.UOO 

590.971 

• 221 

811. 

80A . 

6.75A 

.821 

2 

ACA.UUO 

390.971 

.221 

811. 

60A • 

6.75A 

.821 

3 

AOA.UOO 

390.971 

• 221 

811. 

eoA . 

6.759 

.821 


OUTFLOW EOUALS 20.259 RG/HR 


Ft 

P31 

PA 

HA 

TAT 

TA / 

P5T 

P5 

N5 

I5T T5 /TCTIF 

WOUT 

KT 

RT 

RT 

RHOV 

rhovso 

ITRS 

HUW 

tN/CN**2 1 


CK > 

/1N/CN««2I 



(H> / (K> 

(KG/HRl 



CoRR 

RATIO 

RATIO 


1 

390.971 

37A.327 

• 191 

99A . 

987./ 

382. 9A8 

373.AO0 

♦ 195 

1029. 1022. /1U23. 

A .83 3 

• BAD 

.939 

1.000 

1 .26 3 

.768 

2 

2 

390.971 

371.911 

.205 

982. 

975./ 

381.768 

370.600 

.21U 

1020. 1012./101A. 

5.176 

.839 

.937 

1 »ouu 

1 .256 

.730 

2 

3 

390 .971 

359.816 

.217 

971. 

963./ 

360. 726 

366.500 

.222 

1012. 1005./1005. 

5 .967 

• 5 38 

.936 

1.000 

1 . 2A3 

.709 

2 

A 

39U.971 

366. A22 

.2 35 

1001 • 

991 •/ 

379.010 

36A.700 

• 2AZ 

1U9A. 1033./1037. 

9.783 

• B36 

.935 

1.000 

1 .20A 

.687 

2 


HEAT transfer RESULTS 


FC 

HE • 1 

>TRANSFER-C0EF FK lENTS 

H'NOO 'FACTORS 

COOLED 

GAS 

WALL' 

IE9PERA TURE 

AVG.'THERN.'CONO. 

ETA 

ll‘ 

HOW 

H6U 

HGl 

FC-HOLE 

inpG 

FC'HOLE 

iNPG 

AREA 

TENP 

OUTSIDE 

iniface 

INSIDE 

METAL 

COATING 

ITCO-TCI/ 




«J/(HWA2ASEC*K » 1 



(CNA«2 1 

CK • 

(K 1 

txi 

tRI 

(J/ (CMASECAK 1 » 

(TWO-TCI 


1 

5277. 

3972. 

9733 . 

8890. 

I .000 1 

1 .OUO 

• 968 

255U. 

153A« 

1232. 

1132. 

.388 

.019 

.3019 

3 

2 

5816. 

9 256. 

lOIAU. 

8892. 

1.000 1 

1 .000 

.966 

2550. 

I5A3* 

123A . 

1133. 

• 389 

.019 

.2856 

3 

3 

6 38 A. 

AA83. 

10397. 

8882. 

1.000 1 

1.000 

.968 

2550. 

1592* 

1229. 

1129. 

.387 

.019 

• 27A3 

3 

A 

695 1 . 

A767. 

10915. 

8908. 

1.000 ] 

l.OOO 

.968 

2550. 

1562* 

1299 . 

1191. 

.391 

.019 

.3101 

S 


81 



TABLE Vn. - EXAMPLE 2 (BLADE) CHAMBER OUTPUT 


-output FOR CHAMBER 2 


SI StSItH Of UMltS 

COOlANt GAS constant: 287. OSO J/(K6-KI 

THIS CASE INCLUDES CENTRIFUGAL EFFECTS. ROTATIONAL SPEED EQUALS 16B2S.OU RPN. 
THIS CASE IS FLOW analysis ONLY AND INCLUDES NO NEIAL TERPERA1URL CALCULATIONS 


IS ROWS OF IMPINGENENT HOLES 


ROW HOLES DIAKETER IRHI WALL 

THICKNESS 


L/D 


HOLE IttPlNSENENT R1 PIT 

SPACING DISTANCE IRNt tN/CR**2> 


1 2 

2 2 

3 2 

« 2 

5 2 

b 2 

7 2 

H 2 

V 2 

lU 2 

11 2 

12 2 

li 2 

lA 2 

IS 2 


.Ail8 
.A318 
.A318 
.A318 
.A318 
.4318 
.4316 
.431 B 
.4318 
.4318 
.4318 
.4318 
.4318 
.4318 
.4318 


.381 

.381 

.381 

.361 

.381 

.361 

.381 

.381 

.381 

.381 

.381 

.381 

.381 

.381 

.381 


.882 1.81U 
.882 3.810 
.882 3.810 
.882 3.810 
.882 3.61U 
.882 3.810 
.882 3.810 
.882 3.810 
.882 3.810 
.882 3.810 
.882 3.810 
.882 3.810 
.882 3.810 
.882 3.810 
.882 3.810 


.7b2 217.200 284. 3UU 
.762 219.700 286.400 
.762 222.300 288.500 
.762 224.800 290.700 
.762 227.300 293. UOO 
.762 229.900 295. 2UU 
.762 232.400 297.600 
.762 235.000 299. ago 
*762 237. SDO 302.300 
.762 240. UUO 304 . BOU 
.762 242.600 307.300 
.762 245.100 309. SOU 
. 762 247.700 312. 4UU 
.762 250.200 315.100 
•762 252.700 317.800 


IS NOWS OF FILM COOLING HOLES 


ROW 

HOLES 

OlANETCR INHI 

THICKNESS 
WALL COATING 

L/0 

TTOTALI 

MOLE 

SPACING 

ALPHA 

«OCG) 

beta 

tOEG) 

RHOVG 

IKG/M*»2*HRI 

RH0V2G 

IKG/N«HR*«2> 

R4 

INNA 

P6 

IN/CN442A 

1 

2 

.4572 

1.016 

.000 

4.444 

2.S4U 

iU.UOU 

.OUO 

.50890*07 

.58730*13 

217.200 

264.500 

2 

2 

,4572 

1.016 

.000 

4 .N44 

2.54U 

3U.Q0U 

.000 

.51080*07 

.59020*13 

219.700 

265.200 

3 

2 

.4572 

1.016 

.300 

4.444 

2.540 

30.000 

.900 

.51270*07 

.59300*13 

222 .30U 

266. J9'J 

« 

2 

.4572 

1.016 

.300 

4.444 

2.540 

30.900 

.309 

.51450*0T 

.59580*13 

224. BOO 

266.900 

S 

2 

.4572 

1.016 

.000 

4.444 

2*540 

30.900 

.000 

.51630*07 

.59870*13 

227.500 

267.600 

6 

2 

.4572 

1.016 

.000 

4.444 

2*540 

30.0UU 

.300 

.51820*07 

.60150*13 

229.900 

26B.400 

7 

2 

.4572 

1.U16 

.000 

4.444 

2*540 

3U.UUU 

.OUO 

.52000*07 

,60440*13 

232.400 

269.200 

8 

2 

.4572 

1.016 

• OOO 

4 .444 

2.S4U 

30.00U 

.300 

.52190*07 

.60720*15 

235.000 

269. VUO 

V 

2 

.4572 

1.016 

.000 

4.444 

2.540 

30.000 

.000 

.52370*07 

.61010*1 3 

237.500 

Z7O.700 

lU 

2 

.4572 

1.016 

.000 

4.444 

2*540 

30.000 

.000 

.S?&60«07 

.61290*13 

240.000 

271.500 

11 

2 

. 45 72 

l.OU 

.000 

4 .444 

2*540 

30.300 

.000 

.52750*07 

.61560*13 

242.600 

272.300 

12 

2 

.4572 

1.016 

• 000 

4 .444 

2.54U 

3U.00U 

.000 

.52930*07 

.61860*13 

245.100 

273. 1 JO 

13 

2 

.4572 

1.016 

.000 

4.444 

2.540 

30.000 

.000 

.53110*07 

.62150*1 3 

247.700 

2T3.I00 

14 

2 

.4572 

1.016 

«ouo 

4.444 

2*540 

30*000 

.000 

.53300*07 

.62430*13 

250.200 

274.600 

IS 

2 

.4572 

1.016 

.ouo 

4 ,444 

2*540 

30.0UU 

.000 

.53480*07 

.62720*13 

252. 7UU 

275. 4 JJ 


IHPINGERENt AND FILM COOLING FLOWS HAVE CONVERGED iN 9 OVERALL ITERATIONS 
INFLOW EQUALS 23.452 KG/HR 


IMP 

ROW 

pspltt 

(N/CN**2I 

P2 

H2 

T2T 

(Kl 

T2 

VIHP 

TKG/HRI 

COIHP 


264.300 

271 .933 

.257 

811. 

802 . 

1 .481 

.828 

2 

266.400 

273.921 

.257 

811. 

802. 

1.493 

.826 


2 68*^«0 

276.026 

.256 

611. 

802. 

1 .498 

.828 


290.700 

278.092 

.257 

811. 

802. 

1.512 

*828 

5 

293. UUO 

280.196 

.258 

811. 

802. 

1.530 

.828 

6 

295.200 

282.426 

.256 

811. 

602. 

1 .533 

.828 


297. 6UD 

284 .610 

.257 

811. 

802. 

1 .553 

.828 


299.9U0 

266.926 

.256 

611. 

802 . 

J .558 

.826 

9 

302.300 

289.194 

.257 

611. 

802. 

1.572 

.826 

10 

304. 8U0 

291 .SU5 

.257 

611. 

802 . 

1 .590 

.828 

11 

307. 3UU 

293.954 

.257 

811. 

802. 

1 .599 

.828 

12 

3 uv.euo 

296.354 

.257 

811. 

802. 

1 .612 

.828 

13 

312. 4U0 

298 .896 

.256 

811. 

802 . 

1.622 

.820 

14 

3 15.1uO 

301 .387 

.257 

811. 

802. 

1.642 

.828 

15 

317. BOO 

3U3.924 

.256 

811. 

8 02 . 

1.658 

.828 



OUTFLOW EQUALS 

23.441 

R6/HR 













- 

FC 

P3I 

P4 

N4 

T4T 

T4 / 

P5T 

P5 

N5 

T5T 

T5 /TCTIF 

WOUT 

K 

RT 

RT 

RHOV 

RHOVSO 

IT7S 

RUW 

<N/CH«*2A 


TKI 

/IN/CN**2A 



IKI 

/ 

<RI 

(KG/HRI 



CORR 

RATIO 

RATIO 



271 ,933 

264.642 

. 15U 

411 . 

806./ 

268.533 

264.500 

.150 

811. 

808./ 

0. 

.957 

.843 

.846 

1 .uuu 

.677 

.177 

2 

2 

273 .921 

265. 36U 

.162 

811. 

807./ 

269.933 

265.200 

.162 

811. 

8C7./ 

0. 

I.ObU 

.843 

.864 

l.OOQ 

. 732 

.207 

2 

3 

2 76.026 

266. 178 

.173 

811 . 

607./ 

27] .444 

266. UOO 

.173 

811. 

3U7,/ 

0. 

1.156 

.642 

.877 

l.UUU 

.783 

.237 

2 

4 

278 .092 

266.993 

. 164 

811. 

6U6./ 

272.933 

26 6. SOU 

.184 

811. 

896./ 

9. 

1.242 

.841 

.666 

1.000 

.B3U 

.265 

2 

S 

28U.1V6 

267 .606 

.194 

311. 

BOb./ 

274.444 

267.600 

. 194 

811. 

B06./ 

0. 

1*325 

. 840 

.893 

1.000 

.875 

.295 

2 

t> 

282.426 

266.624 

. 2U4 

811. 

6U5./ 

276. U24 

268. «U0 

.204 

Bll. 

805./ 

0. 

1.411 

.840 

.899 

l.OUU 

,922 

.527 

2 

7 

284 .61U 

269.437 

.214 

81 1 . 

805./ 

277.581 

269.209 

.214 

011 . 

BU4 ./ 

0. 

1.489 

.8 39 

.9U4 

1.000 

.965 

.357 

2 

e 

266 .926 

27U.151 

.2 24 

311 . 

8U4./ 

274.164 

269.900 

.224 

811. 

804 ./ 

0. 

1.577 

.838 

.9U9 

1 .uuu 

1 .U12 

.393 

2 

V 

289 . 194 

270.963 

.233 

611. 

803,/ 

280.769 

270.700 

.234 

811. 

803./ 

9. 

1.65 3 

.837 

.913 

1.000 

1 .U5 3 

.425 

2 

lu 

291 .5U5 

271.774 

. 242 

eii. 

8U3./ 

287.397 

271.500 

.243 

811 . 

803./ 

0. 

1.73U 

.8 36 

.917 

l.UUU 

1 .U94 

.457 

2 

11 

293 .954 

272.592 

.252 

611. 

8U2./ 

284. 103 

272. 30J 

.252 

811. 

602./ 

0. 

1.811 

.835 

.921 

l.UUO 

1 . 136 

.493 

2 

12 

296 .354 

273.401 

.2bU 

Bll . 

801./ 

285. 763 

273.109 

.261 

011. 

801 ./ 

g. 

1.887 

.834 

.924 

1 .000 

1 . 175 

.527 

2 

13 

296 . 696 

274 .111 

.2 7U 

611 . 

601 . / 

287,498 

273. 80U 

.271 

011. 

801 ./ 

9. 

1.972 

.832 

.927 

l.UUU 

1 . 219 

.566 

2 

14 

3UI .367 

274.916 

.2 79 

611 . 

BOO./ 

289.231 

274.600 

.279 

811 . 

800./ 

0. 

2.047 

.8 31 

.9 30 

1.000 

1.258 

.602 

2 

13 

3U3 .924 

275.724 

.267 

311 . 

799./ 

29U.991 

275.400 

.286 

Bll . 

799,/ 

u. 

2. 173 

.8 39 

.933 

l.OUU 

1 .296 

.638 

2 
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Figure 3. - Definitions of film-cooling-hoie angies. 







I- Chordwise film- 
cooling row 


Figure 4. - Vane chamber division. 



cooling rows 

Figure 5. - Blade chamber division. 


84 




Calculate impingement h byeq. (Bill 
(fluid properties dependent on jl 
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Figure 6. - Flow diagram for iterative heat-transfer calcuiations. (Equa- 
tions for coated shell are marked with an asterisk. ) 










Impingement /" 3 
rows; I r 2 



Impinge- 

Number of 

Hole 

ment 

holes per 

diameter, 

row 

row 

cm 

1 

15 

0.0508 

2 

15 

.0508 

3 

15 

.0508 


Film- 

cooling 

row 

Number of 
holes per 
row 

Hole diameter, 
cm 

Ma in- 
stream 
static 
pressure, 

N/cm‘ 

Main-stream 

density 

times 

velocity. 

pv. 

kg/tm^- hr) 

Main -stream 
density 
times 
velocity 
squared, 
PV2, , 
kg/(m • hr'^) 

Main-stream 

HGO,® 

Jf(m‘ ■ sec • K) 

Main-stream 

HG1.» 

J/tm'^ • sec ■ K) 

1 

15 

0.04064 

373.4 

4.364x10^ 

3. 204x10^^ 

5277 

3972 

2 



.04064 

370.8 

4.781 

3.872 

5816 

4256 

3 



.04064 

368.5 

5.107 

4.439 „ 

6384 

4483 

4 



.0381 

364.7 

5.590 

5. 362x10“ 

6951 

4767 


®Main-stream-gas heat-transfer coefficient for coolant temperature equal to main-stream gas temperature. 
°Main-stream-gas heat-transfer coefficient for coolant temperature equal to shell outer temperature. 


Figure 8. - Vane chamber of example 1. 
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Main-stream 
density 
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velocity 
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87 


Figure 9. - Blade chamber of example 2. 
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